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Abstract 



The half-life for two-neutrino double beta decay of Nd has been measured with 
data taken by the NEMO 3 experiment at the Modane Underground Laboratory. 
Using 924.7 days of data recorded with 36.55 g of ^^°Nd the half-life of this process 
is measured to be 

= (9.11+°:! (stat.) ± 0.62 (syst.)) x 10^^ y. 

This result is the world's most accurate half-life measurement of this isotope. A 
search for neutrinoless double beta decay of this isotope is presented using the same 
data taking period. No significant excess of events above the background expectation 
is observed and the limit on the half-life of this process is set to be 

Tf/s > 1.8 X 10^2 y (90% CL). 

This result has significantly improved the previous limit on neutrinoless double beta 
decay of this isotope. Limits are also set on several other neutrinoless double decay 
modes. 
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Chapter 1 
Introduction 



Neutrinoless double beta decay is a lepton number violating process which provides 
unique evidence that neutrinos are Majorana particles, i.e. their own anti-particle. 
Neutrinoless double beta decay is also sensitive to the neutrino mass scale. The 
experimental signature of 0^(3(3 decay is the observation of two electrons, for which 
the total energy sum is equal to the nuclear transition energy. 

This thesis presents a measurement of the half-life of neodymium-150 (^^°Nd) two- 
neutrino double beta decay [2^(3(3) and a search for different modes of neutrinoless 
double beta decay {Qv(3(3). The data used for this thesis were collected by NEMO 3 
between 2003 and 2006, corresponding to 924.7 days of data taking. 

The ^^'^Nd isotope has a nuclear transition energy of 3.367 MeV, which is higher 
than for most of the natural radioactive sources of background, and has a large 
phase space factor. These properties have made ^^°Nd a strong candidate for next 
generation double beta decay experiments, such as SuperNEMO [1] and SNO-I- [2]. 

In NEMO 3 the backgrounds to Qv(3(3 are divided into 2v(3(3 decay and radioactive 
backgrounds. The 2^(3(3 background is irreducible as it has the same event topology 
as Qv(3(3 decay. The precise half-life measurement of 2v(3(3 decay is therefore impor- 
tant for Qvj3j3 searches. It also helps to improve the understanding of nuclear matrix 



9 



elements (NME), which are the major source of uncertainty in the derivation of the 
neutrino mass from the Oz//5/5 half-life. The radioactive backgrounds are reduced by 
applying two-electron event selection criteria. The measurement of their activities is 
necessary for estimating the number of remaining events due to these backgrounds. 
This is achieved by studying control channels with final states different from the 
signal. 

The thesis is set out as follows. The theoretical background to the work presented in 
this thesis and the current status of double beta decay experiments are described in 
Chapter 2. Chapter 4 describes the NEMO 3 experiment. The analysis techniques 
used are discussed in Chapter 5. The radioactive backgrounds to double beta decay 
of ^^°Nd are introduced and their activity measurements are described in Chapter 6. 
This chapter also shows that the background estimation can describe data well in 
several different analysis channels. 

Chapter 7 is dedicated to a measurement of the ^^°Nd 2vf3(3 half-life. The systematic 
uncertainty on the measurement is estimated and the result is compared to other 
measurements of the ^^''Nd half-life. Chapter 8 presents the limits on half-lives 
of different neutrinoless double beta modes and compares these results with other 
searches for new physics in double beta decay experiments. Finally, Chapter 9 
summarises the work described and gives a conclusion. 
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Chapter 2 



Theoretical background 

2.1 Introduction 

The neutrino was first proposed by Wolfgang Pauli in 1930 [3] as a liglit neutral 
particle to solve tlie observed non-conservation of energy in beta decay. Since then 
much work has been done to estabhsh a theory that can describe the fundamental 
particles and the observed interactions between them. As a result, neutrinos have 
become one of the building blocks of the Standard Model (SM) of particle physics, 
but several of their properties such as mass and their Dirac or Majorana nature are 
not known. 

This chapter presents the theoretical background for the measurements and searches 
presented in this thesis. It begins with a brief introduction to the Standard Model. 
Sections 2.3 to 2.4 review the properties of neutrinos with emphasis on the properties 
which are not known. These sections also explain how 01/(3(3 decay can answer some 
of the questions regarding the nature and masses of the neutrinos. Section 2.5 
gives details of the different double beta decay theories and modes. Section 2.5.5 
introduces the nuclear matrix elemenet of the double beta decay. 



11 



2.2 The Standard Model 



In the SM there are two general classes of fundamental particles: fermions, which 
have non-integer spin; and bosons, which have integer spin. The twelve types of 
fermions are subdivided into two groups, leptons and quarks. Leptons have three 
flavours and consist of the charged electron, muon and tau, together with three 
corresponding charge-neutral neutrinos: electron neutrino, muon neutrino and tau 
neutrino. The interactions of particles in the SM are mediated by the exchange 
of gauge bosons. There are three types of interaction in the SM: electromagnetic, 
weak and strong interactions. The properties of fermions and gauge bosons in the 
SM are given in Tables 2.1 and Table 2.2, respectively. In the SM fermions and 
gauge bosons obtain masses through the Higgs mechanism [4] which introduces an 
additional field with an associated particle, the Higgs boson. 

Definition of lepton number in SM 

Each generation of leptons has associated with it a quantum number. There are 
three lepton numbers: electron number (Lg), muon number (L^) and tau number 
(Lr). The electron and electron neutrino have Lg = 1 and the positron and electron 
anti-neutrino have = —1. For all other leptons Lg = 0. Similarly, = and 
Lt- = for leptons not in the muon and tau generations, respectively. In the SM, 
the sums of lepton numbers are conserved in all known interactions. 

2.3 The nature of massive neutrinos 

Despite the success of the SM, it is not able to account for massive neutrinos. In the 
SM, neutrinos are considered to be left-handed and this precludes the possibility that 
neutrinos have mass, as helicity is conserved only for a massless particle. However, 
the observation of the oscillation between the different flavours of neutrinos [6, 7] has 
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fermions 


charge 


mass 


leptons 


electron (e) 


-1 


0.51 MeV 




electron-neutrino [ug) 





< 2 eV 




muon (/i) 


-1 


105.6 MeV 




muon- neutrino (z/^) 





< 2 eV 




tau (r) 


-1 


1777 MeV 




tau-neutrino (ur) 





< 2 eV 


quarks 


up (u) 


+2/3 


1.5-4 MeV 




down (d) 


-1/3 


4-8 MeV 




charm (c) 


+2/3 


1.15-1.35 GeV 




strange (s) 


-1/3 


80-130 MeV 




top (t) 


+2/3 


174 ± 5 GeV 




bottom (b) 


-1/3 


4.1-4.4 GeV 



Table 2.1: The three generations of fermions in the SM, their electric charge and 
mass [5]. 



gauge bosons 


interaction 


mass (GeV) 


7 (photon) 


electromagnetic 





Z boson 


weak 


91.188+0.002 


boson 


weak 


80.425 + 0.038 


g (gluon) 


strong 






Table 2.2: The gauge bosons of the SM and their masses [5]. 

shown that the neutrino mass eigenstates (z/i, 1^2 and u^) and the flavour eigenstates 
(i^e; '^fj. and Ut-) are different, and thus neutrinos have mass. This leads to two 
fundamental questions: what is the nature of the neutrino mass and what is the 
mass scale of the neutrino? The latter is discussed in Section 2.4.2. This section 
focuses on the question of the nature of the neutrino. Neutrinos are either Dirac 
particles, which means that they are distinct from their own anti-particles, or they 
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are Majorana particles, meaning they are their own anti-particles. 



2.3.1 Dirac neutrino 



The Dirac mass term of the neutrino Lagrangian is written as [8] 



(2.1) 



where m/) is the non-diagonal Dirac mass matrix, and z/^ and i^^ are the chirally 
right-handed and left-handed flavour neutrino fields defined as 
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eR 



V, 



liR. 



(2.2) 



The second term in Equation 2.1 is the Hermitian conjugate (h.c.) of the first term. 
As right-handed neutrinos do not couple through the weak interaction, the Dirac 
model predicts three right-handed sterile neutrinos which do not interact and can 
not be detected. Dirac mass terms require four independent components 
l/l, z/^j) and thus suggest that neutrinos and anti- neutrinos are distinct and lepton 
number is conserved. The Dirac model leads to small coupling of the neutrino to 
the Higgs field (in comparison with the coupling of other leptons). The Dirac model 
can not explain why the neutrino mass is much lower than that of other leptons. 



2.3.2 Majorana neutrino 

The Majorana mass term may be constructed out of solely the left-handed or right- 
handed neutrino. However, the left-handed mass is not invariant under the elec- 
troweak gauge group. Thus, it is more common to see the mass term in the right- 
handed construction. The neutrino mass Lagrangian can be written as [8] 

C = ^mnM'Kri^R + h.c., (2.3) 
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where m^iM is the right-handed symmetric Majorana mass matrix, and (t'^)^ is the 
charge-conjugate field of z/^, which satisfies the Majorana condition {upiy = Cu'^'^, 
where C is the charge conjugation matrix, and T denotes transposition. This field 
also has the property that {i^^Y = Pl{i^^Y, where Pl is the left-handed projection 
operator, that is {i^^Y is left-handed. This can be considered in terms of the picture 
of moving into a reference frame where the helicity of the massive neutrino fiips. 

The Majorana neutrino field is self-conjugate and therefore Majorana neutrinos are 
their own antiparticles. Among all fermions only neutrinos can be Majorana particles 
as they do not carry charge. Since the Majorana field has only two independent 
components (z/^, u^) compared to the four in the Dirac field, the Majorana theory 
is simpler and more natural than the Dirac theory. As there is no reason for the 
Majorana term to be zero, neutrinos are assumed to be Majorana particles in most 
theories beyond the SM. 

The see-saw mechanism 

The see-saw mechanism [9] is a Dirac-Majorana model which is developed in order to 
justify the smallness of the mass of the neutrinos in comparison with other fermions 
in the SM. If the neutrino has a Dirac mass as other fermions in the SM do and 
also a right-handed Majorana mass defined through Equation 2.3, then its total 
Dirac-Majorana neutrino mass term Cmass is constructed as 




(2.4) 



y niD rriRM ) 

tion 2.4 can be diagonalised to give the mass eigenstates with eigenvalues 



The matrix M. 



V 




is referred to as the neutrino mass matrix. Equa- 




(2.5) 
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If it is assumed that the Dirac mass, m^) of the neutrino is of the same order of 
magnitude as the Dirac mass of other fermions in the SM and that the right-handed 
Majorana mass term is much heavier than the Dirac mass term (m^jv/ ^ f^o), the 

2 

first approximation of these eigenvalues gives Mi = and M2 = mRM- Thus 
this model predicts two physical neutrinos: the light left-handed neutrino and a very 
heavy right-handed Majorana neutrino. If the heavy neutrino mass is at the GUT 
(Grand Unified Theory) scale (~ 10^^ GeV), then the mass of the light neutrino is 
order of few eV. Therefore this model can explain why the mass of the SM neutrino 
is much lower than the mass of the other fermions. As the Majorana mass term 
appears in both neutrino mass eigenvalues, this mechanism predicts that neutrinos 
are Majorana particles. 



2.4 PMNS matrix 



In the full three neutrino mixing framework the weak eigenstates can be expressed as 
superpositions of three neutrino mass eigenstates (mi, m2, m^) linked via a unitary 
matrix U: 

/ ^ / f/el Ue2 Ue3 

U^l U^2 U^3 
\ Url Ur2 Ur3 ) 



This can be rewritten as 



\ / .i\ 



V2 
V ^3 / 



(2.6) 



(2.7) 



This unitary mixing matrix is known as the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
matrix [10, 11]. The unitary matrix U can be parametrised in the following form 



U 



C12C13 



"S12C13 



-S12C23 - Ci2S23Si3e*'^ C12C23 - Si2S23Sl3e*'' 



S23C13 



(2i 



\ ■Sl2'523 " Ci2S23Si3e*'' " C12S23 — 512023^136*'^ C23C13 J 
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D is the diagonal matrix defined as 



/ 1 \ 



D 



e"^! 



(2.9) 



\ e'^^ J 



where Sij, Cij are the sine and cosine of the mixing angle 6ij, 6 is the Dirac phase and 
01 and 02 are the Majorana phases, which only affect Majorana neutrinos. Thus 
the matrix D only appears in the PMNS matrix if neutrinos are Majorana particles. 

Dirac and Majorana phases 

The Dirac and Majorana phases in the PMNS matrix are CP-violating phases. The 
Dirac CP-violating phase, 6, can be measured by oscillation experiments via compar- 
ing the probability of neutrino and anti-neutrino oscillations from a certain flavour 
to another. However, the oscillation experiments are not able to measure the Ma- 
jorana CP-violation phases (pi and 02 as they only appear on the leading diagonal 
of the PMNS matrix and therefore cancel in all measurable quantities in oscillation 
experiments. Neutrinoless double beta decay (Oz//5/?) experiments may be able to 
provide a possible constraint on CP- violation phases, as Ouf^f] is sensitive to the 
value of the effective neutrino mass: 



(m,) = I J2 Uliml = 144^1 + sl,cl,m2e-'"^' + sl,mse-'''^^'-'^\, (2.10) 



where mj = z/j. More details about neutrinoless double beta decay are given in 
Section 2.5.3. 

2.4.1 Neutrino mass hierarchy problem 

The probability for a neutrino to change from one flavour to another is related to 



Based on the observations from neutrino oscillations, various mass models have been 



the difference between the masses squared of the mass eigenstates, Am^, 
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Normal 



Inverted 



> 



C/3 



m, 



-2.4x10" 



-7.6x10" 



m, 



^3 



-7.6x10" 



-2.4x10" 



Figure 2.1: A schematic view of the hierarchy models consistent with data from 
oscillation experiments. The question mark illustrates that the absolute mass scale 
of neutrinos is not known. 

proposed [12]. These can be categorised as the normal hierarchy in which mi has 
the lowest mass among the three, and the inverted hierarchy in which has the 
lowest mass. In both scenarios mi and m2 have similar masses. The third model is 
degenerate, where the three mass eigenstates have similar masses. Figure 2.1 shows a 
schematic view of the normal and inverted hierarchy models. The approximate mass 
squared difference between the mass eigenstates found by oscillation experiments 
is shown [13, 14]. Figure 2.2 shows the distribution of {m^) with respect to the 
oscillation parameters for normal (red) and inverted hierarchy (green). The next 
generation double beta decay experiments will be sensitive to the full degenerate and 
inverted hierarchy regions. If no neutrinoless double beta decay event is observed 
by these experiments, then the inverted hierarchy can be excluded. 



2.4.2 Absolute mass scale of neutrino 

Oscillation experiments are trying to measure the difference between mass eigen- 
states of the neutrinos, but they are not sensitive to the absolute mass scale of the 
neutrino. Experimentally there are three possible ways to constrain or measure the 

18 




lightest neutrino mass in eV 



Figure 2.2: Bounds on the effective neutrino mass with respect to the hghtest neu- 
trino mass at 90% confidence level [12]. The normal hierarchy is shown in red and 
the inverted hierarchy is shown in green. The degenerate region is where the green 
overlaps with the red. The darker regions of the plot show how the bounds on the 
effective neutrino mass would be constrained if the current oscillation parameters 
were calculated with negligible error. Also shown in grey are the regions disfavoured 
by current Oi//3/3-decay limits and from cosmology. 

absolute mass scale. One way is to search for the neutrino rest mass by explor- 
ing the endpoint energy of the electron spectrum in tritium beta decay and thus 
find the mass scale of the neutrino {{rrii,)) directly from the kinematics of the de- 
cay. Currently a limit for the electron neutrino mass of less than 2.2 eV has been 
achieved [15, 16]. The advantage of this experimental approach is that it is indepen- 
dent of the nature of the neutrino (Dirac or Majorana); however, the measurement 
of the mass from the decay spectrum is limited by the energy resolution of the ex- 
periments and background considerations. This makes it experimentally difficult to 
improve the sensitivity of the experiment to an energy scale of meV. 
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The second way to constrain the neutrino mass is from analysing cosmological data. 
The astrophysical experiments are able to set limits on the sum of the neutrino 
masses (X^i"^*)- These limits are based on the cosmological models which give the 
neutrino contribution to the density of the universe. The most recent limits for 
rrii range from 0.7 to 2 eV at 90% confidence level [17]. This range can constrain 
the lightest neutrino mass (mi or m^). The disfavoured region of lightest neutrino 
mass is shown by a vertical gray bound in Figure 2.2. As the cosmological models 
are based on several assumptions, it would be challenging to improve these results. 

The third approach to constrain the absolute mass scale of neutrinos is searching 
for neutrinoless double beta decay {Ou/S/S). As mentioned previously in Section 2.4, 
neutrinoless double beta decay is sensitive to the effective neutrino mass {m,y) if the 
neutrino is a Majorana particle. 

2.5 Theory of double beta decay 
2.5.1 Beta decay 

Beta (/?) decay is a type of radioactive decay in which the weak interaction converts 
a down quark of the neutron into an up quark of a proton while emitting an electron 
and an antineutrino. Thus, this process changes the atomic number, Z, of a nucleus 
by one unit, while the atomic mass A remains the same: 

{A, Z) ^ {A, Z + 1) + e- + (2.11) 

There are also two other possible modes of this process which are known as (3~^ and 
electron capture (EC) and defined as: 

{A,Z) ^ {A,Z -l) + e+ + Ue (^+ decay), (2.12) 
{A,Z) + e- ^ {A,Z -1) + (EC). (2.13) 
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2.5.2 Double beta decay 



Double beta decay, 2u[3[3 [18], is a process in which two beta decays occur simulta- 
neously with emission of two electrons and two antineutrinos: 

{A,Z) = {A,Z + 2) + 2e- + 2u. (2.14) 

Double beta decay can only occur in nuclei with even atomic number and atomic 
mass (even-even). In these nuclei the single beta decay can be either energetically 
forbidden or strongly suppressed. The nuclear transition energy, Qpp, for this process 
is defined as: 

Q0P = m{A,Z) -m{A,Z + 2) -2me, (2.15) 

where m{A, Z) and m^A, Z + 2) are, respectively, the masses of the initial and 
final nuclei and me is the mass of the electron. Figure 2.3 shows a diagram of the 
energy requirements for double beta decay. This process is the rarest known kind of 
radioactive decay, and it is predicted in only 36 isotopes. Double beta decay occurs 
within the SM and involves a second order weak interaction as shown in Figure 2.4a. 
The rate of the process is characterised by its very long lifetime (more than 10^^ 
years) . 

The half-life {Ty2) '^'^ double beta decay process is expressed by the equation: 

(Ti%)-i = G''^(Q,Z)|M2'^|2, (2.16) 

where G'^'^lQpf^, Z) is an analytically calculable phase space integral which scales 
with Q^Jj and the atomic number, Z [19]. The term M'^'^ is the nuclear matrix 
element (NME) (Section 2.5.5) for double beta decay which gives the probability 
of the process and can be calculated theoretically. To test the methods used to 
calculate NMEs, it is important to determine the half-life of 2z//3/3 experimentally. 
Two neutrino double beta decay is an irreducible background for searches for beyond 
the Standard Model (BSM) double beta decay processes for detectors with weaker 
energy resolutions. 
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Forbidden 



^ Z 



Figure 2.3: A schematic diagram of the energy requirements for double beta decay. 
A parent nucleus decays to a daughter nucleus because the intermediate 

single beta decay to z+v^ is energetically forbidden. 

2.5.3 Neutrinoless double beta decay 

Neutrinoless double beta decay [OuPP) [20] involves a transition of two neutrons 
into two protons with the emission of two electrons and no neutrinos (Figure 2.4b): 



(A, Z) = {A,Z + 2) + 2e-. 



(2.17) 




(a) 




n 



W 



(b) 




-> e" 



Figure 2.4: Two-neutrino double beta decay and neutrinoless double beta decay 
diagrams. 
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In this process a right-handed antineutrino emitted at one vertex undergoes a hehcity 
flip and is absorbed into the other vertex as a left-handed neutrino. This can only 
happen if neutrinos and antineutrinos are identical (Majorana particles) and have 
mass. This process violates total lepton number by two units and is forbidden in 
the SM. Observation of neutrinoless double beta decay would prove that neutrinos 
are Majorana particles and would answer the most fundamental question about the 
nature of the neutrinos. 

As there are no neutrinos in the final state, the experimental signature of Oz//3/3 is two 
electrons, with energies summing to the nuclear transition energy (Qpis). For 21/(3(3 
decay the energy sum of the two electrons forms a continuous spectrum (Figure 2.5). 

The half-life of Qv(3(3 is expressed by: 



where {my^) is the effective Majorana neutrino mass and me is the electron rest 
mass. Thus the experimentally measured half-life, or in case of non-observation of 
the process its lower limit, can provide the effective neutrino mass or an upper limit 
on it. The phase space integral, G^'^ , is proportional to and the atomic number, 
Z [21]. Thus the neutrinoless double beta decay rate is higher for isotopes with high 
(5/3/3 value. 

Right-handed current 01//3/3 

The 0^(3(3 mechanism that involves the helicity flip of the neutrino due to its mass 
is known as the mass mechanism. This mechanism is not the only possibility that 
leads to emission of two electrons in the final state. There are theories beyond 
the Standard Model which predict the existence of pure right-handed W bosons 
(Wr) or W bosons with mixtures of left-handed and right-handed W bosons {W = 
Wr + Wl) [19]. If Wpt exists and the neutrino is a Majorana particle, then a right- 
handed neutrino can interact at the other vertex without need for a helicity flip. This 




(2.18) 
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Figure 2.5: Energy sum of the two electrons for 2u/3/3, Oz//5/5, Majoron mode one 
(n=l), Majoron mode two (n=2), Majoron mode three (n=3) and Majoron mode 
seven (n=7) for an arbitrary isotope. 

mechanism is known as the right-handed current mechanism of neutrinoless double 
beta decay. In this thesis a limit is set only on the half-life of pure right-handed 
neutrinoless double beta decay. The half-life of this process is inversely related to 
the phase space factor of this process, the NME {M^'^) and A, the coupling constant 
of the right-handed neutrino with Wr, 



2.5.4 Double beta decay with emission of Majorons 

Several beyond the SM models exist in which global B-L (Baryon-Lepton) symmetry 
is broken spontaneously. These lead to the prediction of a massless Goldstone boson, 
called the Majoron [22] which can couple to the neutrino. There are three types of 
Majoron models: singlet Majoron, doublet Majoron and triplet Majoron [23]. The 




(2.19) 
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measurement of the Z boson width at LEP [24] has ruled out the doublet and 
triplet Majoron models, but singlet and dominantly singlet Majoron models are still 
possible. These models predict that the two neutrinos in double beta decay couple 
with a Majoron, ■)^: 

{A,Z)^{A,Z + 2) + 2e + x\ (2.20) 
Figure 2.6a shows a diagram of such an emission. The drawback of the singlet Ma- 
joron models is that in these models the Majoron couples to the neutrino with a 
coupling strength oi g = {m^j^/MB^L) [25] where m^j^ is the mass of the light neu- 
trino and Mb-l is the symmetry-breaking scale (which is higher than the electroweak 
scale ~ 90 GeV). In order to have an observable rate for Majoron emitting double 
beta decay (Oz//3/9x°); the singlet Majoron model requires severe fine tuning. To 
avoid the fine tuning problem, several new Majoron models have been constructed. 
Here the term Majoron means light or massless boson with couplings to neutrinos. 
In this definition the Majoron is not constrained to be a Goldstone boson. 

There is a doublet Majoron model in supersymmetry which predicts neutrinoless 
double beta decay with emission of two Majorons [23]: 

{A,Z)^{A,Z + 2) + 2e + 2x\ (2.21) 

Figure 2.6b shows the neutrinoless double beta decay diagram with the emission of 
two Majorons. In this diagram the Zino, Z, is the fermionic super symmetric partner 
of the Z boson. There are also other Majoron models which predict double beta 
decay with emission of one or two Majorons. Majorons in these models are predicted 
to be vector bosons and carry lepton charge and have mass [26-28]. 

A model for neutrino masses is proposed in the context of large extra dimensions. 
In this model the global B-L symmetry is broken spontaneously by a gauge singlet 
Higgs field in a bulk [29]. This leads to a bulk singlet Majoron which is observable 
in neutrinoless double beta decay. 

There are in total 10 Majoron models (one bulk Majoron and 9 light or massless 
bosonic Majorons) that are of interest in double beta decay experiments. The energy 
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Figure 2.6: Diagram of double beta decay with emission of one Majoron in the singlet 
Majoron model and of two Majorons from the supersymmetric Majoron model. 

of the final state electrons from these ten models can form four distinguishable 
distributions [27]. The possible two electron energy spectra (Esum) for different 
Majoron modes of ^^°Nd are shown in Figure 2.5. In the figure, the index n defines 
the shape of the spectrum by modifying the distribution with a factor {Qpj3 — EsumY ■ 
Majoron mode one {n = 1) denotes the theories which lead to emission of one 
Majoron, Majoron mode two (n = 2) denotes the bulk Majoron emission, Majoron 
mode three (n = 3) denotes the emission of one or two massless lepton number 
carrying Majorons and Majoron mode seven [n = 7) denotes the emission of two 
light Majorons. Note that electrons from the n = 7 mode have on average lower 
energies than 2i^/?/? electrons; thus, it is more difficult to extract experimentally, 
since most backgrounds tend to dominantly produce low-energy electrons. The 
half-life of Majoron-emitting double beta decay is expressed by 



I-' 1/2 ) 



(2.22) 



The index a indicates the Majoron mode, (ga) is the effective coupling constant 
of Majoron mode a, G^^ is the phase space factor and M^^ is the nuclear matrix 
element of each Majoron mode. 
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2.5.5 Nuclear Matrix elements 



If the half-life of neutrinoless double beta decay is measured, the effective neutrino 
mass can be calculated from Equation 2.18. Accurate calculations of the NME are 
necessary in order to measure the effective neutrino mass with minimum uncertainty. 

Two main theoretical methods have been widely used to calculate the nuclear matrix 
elements for neutrinoless double beta decay: the nuclear shell model (NSM) and 
the quasi-particle random phase approximation (QRPA). The NSM [30] is useful 
for calculating single particle states close to the Fermi level and is only reliable for 
light nuclei such as ^^Ca, ^^Ge and ^^Se. For medium to heavy double beta decay 
isotopes the shell model calculations are difficult to carry out. 

For heavier nuclei the proton- neutron QRPA (pn-QRPA) [31] has been found to be 
a powerful model. The QRPA can handle a great number of intermediate states. In 
double beta decay the initial nucleus decays to the final nucleus through virtual exci- 
tations of all states of the intermediate nucleus. The intermediate excited states are 
obtained by solving the QRPA equations. The QRPA calculations contain two two- 
body interaction matrix elements: particle-particle (pp) which is correlated to the 
proton-proton interaction and particle-hole (ph) which is correlated to the proton- 
neutron interaction. Both matrices contain independent interaction constants Qph 
and gpp [32]. 

The ph matrix elements mainly affect the giant Gamow- Teller resonance [33], which 
is reproduced accurately by the QRPA and as such gph is generally fixed. The gpp 
constant, which has a large effect on NME and double beta decay, is a free parameter. 
This parameter causes uncertainties in NME calculations and thus much of the work 
into QRPA has been focused on the gpp problem. One method in the early stage of 
development is to fix the value of gpp by fitting it to the available 2z//5/5 and (3 decay 
data. Another important issue in QRPA is that the nuclei of experimental interest 
are assumed to have spherical symmetry. As many heavy nuclei are deformed [34] 
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isotope 




™Ge 


3.33-4.58 


82Se 


2.01-4.17 


96 2r 


1.01-1.28 


looMo 


2.22 -3.53 


130 rp J 


2.27-3.77 


150Nd 


3.14-4.04 



Table 2.3: Oz//3/3 NME (M°^) for several double beta decay isotopes, using the QRPA 
method of calculation [35]. 

this approximation does not give a realistic calculation for these nuclei (such as 
^^°Nd). Table 2.3 presents the most recent calculated Qv(3(3 NMEs [M^'^) for several 
double beta decay isotopes. 
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Chapter 3 

Double beta decay experiments 



The challenging task in double beta decay experiments is to search for the Oz//5/3 
peak in the electron energy spectrum within the continuum of 2i^/5/5 and natural 
radioactive background. This chapter presents the experimental criteria and status 
of double beta decay experiments. 



3.1 Experimental criteria 

In order to search for OujSf] signal and minimise backgrounds several factors need to 
be considered in Oz//3/3 experiments [36]: 

• The double beta decay isotopes and the detector components must have high 
purity. This goal can be achieved via purification of double beta decay sources 
and constructing smaller and more granular detector components. 

• To suppress cosmic ray backgrounds the detectors are situated underground. 

• The isotopes with large Q/g/? are favoured as the energy region of interest is 
above potential backgrounds (more details are given in Section 3.1.1). 

• A good energy resolution is required to prevent the tail of the 2z//5/5 spectrum 
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from extending into the Oz//5/5 region of interest. This produces a good signal 
over background ratio. 

• To reject background, event reconstruction and good particle identification is 
required. 

• The NME is understood for some isotopes more than others. The interpreta- 
tion of limits or, in case of discovery, 01/(3(3 signals requires measurements for 
a range of isotopes. 

• Several Oi'f3f3 modes can be studied through energy and angular distributions 
of the electrons in the final states. 

So far, no experiment searching for 01/13(3 has managed to include all the above 
criteria for its 0^(3(3 searches. Section 3.2 gives details of the current experimental 
status of double beta decay detectors. 

3.1.1 Choice of isotopes 

As mentioned in Section 2.5.2, 21/(3(3 occurs in 36 isotopes. However, it is not 
experimentally favourable (and interesting) to detect the decay in all these isotopes. 
Most of these isotopes have a low Q/j/j value which is similar to the Q values of other 
natural radioactive decays. As double beta decay is rare, it will be difficult to detect 
it above background. Currently there are nine double beta decay isotopes which are 
experimentally considered for neutrinoless double beta decay searches: ^'^Ca, ^^Ge, 
^^Se, ^^Zr, i°0Mo, ^^''Cd, ^^oxe, ^^^Xe and ^^^Nd. Table 3.1 gives the Q^^ value and 
natural abundance of each of these isotopes. This thesis is about the double beta 
decay study of ^^°Nd, thus the rest of this section describes the properties of this 
isotope. 
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Transition 


Q/3/3(keV) 


Natural Abundance (%) 


76Ge ^76se 


2039 


7.8 


i36xe ^^^^Ba 


2479 


8.9 


laorpg ^i30xe 


2533 


34.5 


iiecd ^ii6gj^ 


2802 


7.5 


82se ^82Kr 


2995 


9.2 


looMo ^looj^u 


3034 


9.6 


962r ^96Mq 


3350 


2.8 




3367 


5.6 


48Ca ^48rpi 


4271 


0.187 



Table 3.1: Isotopes used for double beta decay studies in experiments. The nuclear 
transition energy (Q/3/3) and the natural abundance of each isotope are shown [19]. 

3.1.2 Neodymium-150 

The double beta decay of ^^°Nd to the ground state of ^^°Sm is accompanied by 
the simultaneous emission of two electrons. The 2uj3l3 decay to the ground state 
of ^^''Sm has already been observed by several experiments. Figure 3.1 shows a 
decay scheme of this isotope. In addition, ^^°Nd can decay to ^^°Sm 0^ and 2^ 
excited states and then de-excite by emitting photons. In the case of a decay to 
the 0^ excited state, the two electrons in the final state are accompanied by two 
photons with energies 410 keV and 334 keV, and in the case of decaying to the 2^ 
state, electrons are emitted with one 334 keV photon. The 2z//?/? half-life of these 
processes is predicted to be greater than 10^° y [37]. A precise measurement of the 
21/(3(3 decay of this isotope to the ground state is important since it is an irreducible 
background to Oiy(3(3 and also can be used to improve current NME calculations. 

The Neodymium-150 (^^'^Nd) isotope has the second highest Q/?/? value of all the 
double beta decay isotopes (after ^^Ca). This, combined with a high atomic num- 
ber, makes the neutrinoless double beta decay phase space factor G^^{Qf^p, Z)) the 
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isoNd 




Figure 3.1: A decay scheme of ^^°Nd to ground state of ^^°Sm. 

highest of all isotopes. Figure 3.2 shows approximate G^'^{Q(^i3, Z) values for several 
double beta decay isotopes. From Equation 2.18 it is shown that the half-life of Qv(3(3 
process is inversely related to G^'^{Q(^i3, Z), thus the event rate for ^^°Nd is expected 
to be higher than all other isotopes, assuming identical NMEs. This feature makes 
^^°Nd a strong candidate for next generation double beta decay experiments such as 
SNO+ [2] and SuperNEMO [1] (Section 3.2). Currently, the only drawback of using 
isoNd for 01/(3(3 searches is the high uncertainty on the NME calculation due to its 
mass and deformation [38]. 




'^Ca '^Ge ''Se ''Zr ""Mo "'Cd ""Te "'Xe '^'NcI 



Figure 3.2: The approximate phase space factor of several double beta decay iso- 
topes [39]. 
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3.2 Experimental status 



The observation of neutrinoless double beta decay can answer important questions 
in neutrino physics. In the past ten years several experiments have been therefore 
constructed and many more are currently in the R&D stage. The double beta detec- 
tors are divided into two main categories: homogeneous (source equal to detector) 
experiments in which the double beta decay source is part of the active detector 
and heterogeneous (source not equal to detector) experiments in which the source 
is separated from the detector. The homogeneous detection technique gives good 
energy resolution. The heterogeneous experiments contain tracking detectors which 
give them a good particle identification. 

3.2.1 Experiments following the homogeneous system 

Semiconductor germanium (Ge) detectors are among the most popular double beta 
decay experiments. This is because of the cheap enrichment process for ^^Ge. The 
NME of this isotope is also relatively well known. There have been two previ- 
ous semiconductor ^^Ge detectors which produced important Ou^f] search results: 
Heidelberg-Moscow and IGEX. Heidelberg-Moscow was a semi-conductor germa- 
nium detector, enriched to 86% in ™Ge. The final design of the detector consisted 
of five individual sub-detectors with a total mass of 11.5 kg. The experiment ran 
between 1990 and 2003. A claim for discovery of Ou/Sf] was made in 2001 by a sub- 
group of the Heidelberg-Moscow collaboration [40]. The half-life value of Oz//5/3 was 
obtained to be = 1-19 ^aso xlO^^ years, with (m^) = 0.2 - 0.6 eV. Figure 3.3 
shows the result of this experiment. However, this result has received criticism [41]. 
It is believed that the background and the systematic uncertainties of this experi- 
ment are underestimated. The IGEX experiment used a similar experimental tech- 
nique to the Heidelberg-Moscow experiment and produced a limit on the half-life of 
T°;2 > 1.57 X 10^5 years at 90% confidence level (CL) [42]. 
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Figure 3.3: Energy sum of the two electrons in Ou/SP region for "^^Ge isotope 
of Heidelberg-Moscow experiment. A peak is observed around the value of 
^•^Ge (2.039 MeV) [40]. 

To confirm or disprove the Heidelberg-Moscow results, two experiments are being 
designed to study the region highlighted by the Heidelberg-Moscow experiment: 
GERDA [43] and MAJORANA [42]. Both of these detectors will use Ge enriched 
to 86% ''^Ge. In the first phase of its running, GERDA will utilize the existing 
Heidelberg-Moscow and IGEX experiments (18 kg of ^^Ge) and will reach a half-life 
sensitivity of 3 x 10^^ years. This sensitivity can rule out or confirm the Heidelberg- 
Moscow results. In the second phase of its running another 38 kg will be added, 
giving a half-life sensitivity of ~ 1.4 x 10^^ years, which corresponds to (mj,) ~ 
124 meV. MAJORANA's eventual half-life sensitivity is hoped to be 5.5 x 10^*^ y, 
corresponding to {rriiy) ~ 61 meV. The final design of this experiment is still under 
development. 

One running source equal detector experiment is CUORICINO [44]. This experiment 
uses bolometers containing double beta decay isotopes, running at extremely low 
temperatures (order of mK). When a double beta decay event occurs, the electrons 
of the decay increase the temperature by depositing energy inside a crystal. By 
measuring the temperature rise the deposited energy is measured. Each bolometer 
is 5 X 5 X 5 cm^ and is fixed into a layered tower structure. The total mass of the 
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crystals is 41 kg, corresponding to 11 kg of ^^'^Te. This detector has set a hmit on the 
half-hfe of neutrinoless double beta decay of ^^^Te to be more than 3.0 x 10^^ years, 
corresponding to {m^) < 0.42 - 0.58 eV [44]. CUORE [45], which is the next 
generation of the CUORICINO experiment, will have 760 kg of Te02 crystals with 
a total of 206 kg ^^'^Te. It will consist of 35 towers which have the same size as the 
CUORICINO tower. It is expected to reach a half-life sensitivity of 2.5 x 10^^ years 
which corresponds to a mass sensitivity of {m,^) ~ 45 — 53 meV, depending on NME. 

The COBRA [46] experiment is in the R&D stage and is going to use CdZnTe 
semiconductor crystals. CdZnTe contains nine double beta isotopes, with ^^"Te 
and ^^^Cd having the highest abundance in the detector. This experiment is also 
searching for OuP^jS^, OujS'^ EC and Oi^EC EC processes by studying ^'^^Cd, ^^Zn, 
^°^Cd and ^^'^Te decays. The finished experiment will hold 64000, 1 cm^ CdZnTe 
detectors. The detector will have 418 kg mass, in which 183 kg will be Cd enriched 
to 90% in ^^^Cd. The half-life sensitivity of this experiment will be greater than 
10^^ years, corresponding to {rriiy) ~ 38 — 96 meV. 

The SNO+ experiment is planning to use most of the SNO experiment infrastructure 
and shielding (SNO was a successful experiment that confirmed neutrino oscillations 
in solar neutrinos), replacing the heavy water with 1000 tones of Nd-loaded liquid 
scintillator. A 1% Nd-loading would correspond to 560 kg of natural Nd [47] or 56 kg 
of ^^°Nd. It is also a possibility that this experiment could load the scintillators with 
enriched ^^''Nd. This experiment is multi-purpose and apart from double beta decay 
studies is going to focus on low energy solar neutrinos as well as geo-neutrinos, 
reactor and supernova neutrinos [47]. 

3.2.2 Experiments following the heterogeneous system 

These type of experiments are divided into two categories: the time projection 
chambers (TPCs), where the double beta decay isotope is a filling gas or liquid, and 
tracking-calorimeter experiments, where the double beta decay source is a solid foil. 
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The EXO [48] experiment is a TPC experiment. Its goal is to use between 1 to 10 
tons of Xe enriched to 80% in ^'^^Xe. The experiment is developing a high resolution 
liquid xenon TPC with good tracking capability. A high voltage cathode will be 
located in the middle of the TPC volume. At each end there will be 114 horizontal 
and vertical anode wires. These wires will be for charge collection, allowing energy 
measurement and particle trajectory reconstruction. The liquid xenon is also an 
efficient scintillator. The use of both collected charge and scintillator light improves 
the energy resolution. The experiment will have two phases; the first phase is called 
EXO-200 (200 kg of ^^exe is used) and is expected to reach a mass sensitivity of 
30 meV. In the second phase the EXO detector is expected to trap and identify the 
^^^Ba ions, the daughter isotope of ^^^Xe, by laser spectroscopy. These ions can be 
identified via atomic spectroscopy by optical pumping with blue and red lasers [49]. 
The radioactive background events can be significantly reduced by Ba tagging. 

NEMO 3 [50] follows the tracker-calorimeter technique and is currently the only run- 
ning experiment of this type. More details about this detector are given in Chapter 4. 
The SuperNEMO detector is a next generation double beta decay experiment which 
is going to use the same experimental approach as NEMO 3 but will improve the 
sensitivity to neutrino mass scale by one order of magnitude {{rrii,) < 45 — 85 meV). 
This experiment will comprise 20 individual modules, each of which has a rectan- 
gular shape with the source foil in the centre surrounded by the tracker and the 
calorimeter. This experiment will use 100 kg-200 kg ^^°Nd and/or ^^Se isotope. 

There are also two other experiments which are planning to use heterogeneous tech- 
niques: MOON [51] and DCBA [52]. MOON will use one tonne of ^°°Mo and hopes 
to achieve an effective neutrino mass sensitivity of {m,y) < 50 meV. DCBA will use 
natural ^^°Nd in its first phase of data taking and then enriched ^^°Nd in the second 
phase. They hope to reach a sensitivity of {m^) < 0.12 eV. Table 3.2 summarises the 
current and future double beta decay experiments, the isotopes they use for Oi^/5/5 
searches, the isotope masses, their expected effective neutrino mass sensitivity and 
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Experiment 


Isotope 


mass (kg) 


Expected (mi,)(meV) 


Time scale 


GERDA Phase 1 


^6Gc 


18 


<250-440 


2010 


GERDA Phase 2 


^6Ge 


56 


<124 


2011 


CUORICINO 


i30Te 


11 


<420-580 


current limit 


GUORE 


i30Te 


206 


<45-53 


2011-2016 


SNO+ 


iso^d 


56 or 560 


under study 


2012 


COBRA 


"«Gd 


164 


<38-96 


- 


EXO-200 


i36xe 


200 


<30 


2009 


NEMO 3 


io"Mo 


7 


< 630-110 


current limit 


SupcrNEMO 


i50Nd/82Se 


100-200 


<45-85 


2012 


MOON 


i""Mo 


1000 


<50 


> 2011 


DCBA phase 1 


natural Nd 


600 


<120 


> 2011 



Table 3.2: Neutrinoless double beta decay experiments, the mass of the isotopes 
they use, their expected sensitivity and their running time scale 

their running time scale. 
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Chapter 4 

The NEMO 3 detector 



The Neutrino Ettore Majorana Observatory 3 (NEMO 3) experiment has been run- 
ning in the Frejus Underground Laboratory (Laboratoire Souterrain de Modane) in 
France since February 2003. The main objective of the experiment is to search for 
evidence of neutrinoless double beta decay in a variety of isotopes. 

The NEMO 3 detector is a heterogeneous system, in which the radioactive sources 
do not make up part of the active detector. Particle tracking and energy mea- 
surements are also performed independently, giving good particle identification and 
understanding of radioactive backgrounds. 

4.1 General description of the NEMO 3 detector 

NEMO 3 is cylindrical in design and is segmented into 20 sectors. The sectors are 
numbered from to 19. Each sector contains tracking chambers and calorimeters in 
order to directly detect decay particles emanating from the radioactive source foils. 
The source foils are positioned in the centre of each sector and divide the sectors 
into two parts: the inner part, which is from the source foil to internal wall of each 
sector; and the outer part, which is from the source foil to external wall of each 
sector. The whole detector is surrounded by a solenoid providing a magnetic field 
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of 25 Gauss, external neutron shielding and an anti-radon tent. 



The NEMO 3 detector can be defined in both right-handed Cartesian and cylindrical 
coordinate systems. In the Cartesian system, the x axis starts from the centre of 
the detector and is along the edge of sector 0. The y axis starts from the centre of 
the detector and is along the edge of sector 5. The z axis starts from the centre 
of the detector and points vertically upward. By performing notational conversion, 
the radius {R = + y^), polar angle (</> = arccos(-|)) and z of the cylindrical 
coordinates are formed. Figures 4.1 and 4.2 show a cutaway view of the NEMO 3 
detector and a view of one of the sectors of the detector, respectively. 

central tower 



Iron shield 



tank containing 
borated water 



wood 




wood 



magnetic coll 



Figure 4.1: A cutaway view of the NEMO 3 detector [50]. The magnetic coil, iron 
7-ray shield, and neutron shielding are shown. 
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Figure 4.2: A sector of NEMO 3 viewed from two different angles [50]. The source 
foil, scintillator blocks and photomultipliers are shown. The tracker cells are located 
between the internal and external walls on both sides of the foil. The tracker wires, 
which run vertically, are not shown for clarity. 

4.2 The NEMO 3 sources 

Unlike all other currently running double beta decay experiments where calorimeters 
serve as both double beta decay source and detector, the sources in NEMO 3 are 
independent from the detector. Inside NEMO 3 there are seven double beta decay 
isotopes, thus permitting study and comparison of results from different isotopes. 
Each sector contains seven foils mounted vertically at a radius of 155 cm from the 
centre within a metallic support frame. This structure allows multiple isotopes to 
be mounted within one sector. Each source foil has a height of 2.5 m, a width of 
65 mm and an area density of 30-60 mg/cm^. 

The choice of source nuclei for NEMO 3 was based on several factors: the double beta 
decay transition energy, Qpis, the phase space factors, G^'^ and G'^'^, corresponding to 
neutrinoless decay and two-neutrino decay; the nuclear matrix elements {M^'^ and 
M^'^); the background in the energy region around the Q/j/j value; and the natural 
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Figure 4.3: The source foil distribution in the 20 sectors of NEMO 3. 
abundance and enrichment possibihties of the isotope. 

Nine isotopes in various quantities were chosen for NEMO 3 and these are placed in 
the order and quantities displayed in Figure 4.3. Double beta decay occurs in seven 
of them: ^^^Cd, ^^Se, ^°°Mo, ^^Zr, ^^^Nd, ^^Ca and ^=^°Te. The two other isotopes, 
Cu and natural Te, have a negligible impurity and are mainly used for measuring the 
background from external sources. The ^'^''Mo and ^^Se foils have the highest mass 
in the experiment, thus much effort has been focused by the NEMO collaboration 
into neutrinoless double beta decay searches of these isotopes [53-55] . 

The neodymium source foil 

The neodymium-150 composite foil (Nd203) was enriched (95.0 ±0.5% of isotope 
150) by the Institute for Nuclear Research of Moscow. It is situated in foil 6 of 
sector 5 (Figure 5.1). The length of the active foil is 234 cm and the width is 
6.5 cm. The total mass of the ^^°Nd foil is 57 g, in which 36.5 ± 0.1 g is ^^°Nd, 
6.458 g is a Mylar support film and the remainder is composed of foil impurities and 
polyvinyl alcohol (PVA). The impurities inside the foil were measured using a high 
purity germanium (HPGe) detector to be ^^^""Pa, ^o^Bi, 154^^^ 1522^^2143-^ 2i4p^^ 
208TI and [50]. 
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Figure 4.4: A view of the sector 5 of the NEMO 3. The ^^''Nd foil is shown with 
^o'^Mo foil on its left and ^^Ca and ^^Zr on its right. 

4.3 The tracking detector 

The tracker consists of 6180 octagonal drift cells operating in Geiger mode [56]. 
Each of the 20 sectors has an internal and external tracking volume, each containing 
drift cells arranged in a 4-2-3 layer configuration as shown in Figure 4.5. This gives 
a total of nine drift cell layers on each side of the source foil to reconstruct the 
particle tracks. 

Each cell has a diameter of 3 cm and a length of 270 cm, and contains a central 
anode wire surrounded by eight ground wires. The layers of drift cells are separated 
by an extra ground wire in order to reduce electrostatic cross-talk. All wires are 
composed of stainless steel and have a diameter of 50 fim. The tracking cells are 
strung between top and bottom walls. At the top and bottom of the wires there 
are 3 cm long copper cathode rings with a diameter of 2.3 cm. Figure 4.6 shows a 
diagram of a drift cell in NEMO 3. 
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anode wire ground wire 




Figure 4.5: A top view of a tracking sector of NEMO 3. The 4-2-3 drift cell layer 
configurations are shown in both the external and internal parts [50]. 

Drift cell operation in NEMO 3 

The tracking system is immersed in a gas mixture of helium (95%), argon (1%), ethyl 
alcohol (4%) and water (1500 ppm). The Geiger mode operating voltage of the cells 




Figure 4.6: A drift cell from NEMO 3 [50]. The anode wire, cathode rings and the 
ground wires are illustrated. 
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is around 1620 V. When a charged particle passes through a cell, it causes ionisation 
of the gas atoms. The result is a track of He"'" ions and free electrons. These electrons 
drift towards the anode wire and cause further ionization, triggering an avalanche 
process as the electrons accelerate. The avalanche process becomes significant in the 
region very close to the wire (around 1 mm), where the electric field is high. When 
the avalanche arrives on the anode, it causes a pulse which has a typical rise time of 
10 ns. The time to digital convertor (TDC) connected to the anode gives the arrival 
time of the initial avalanche to the anode and thus the transverse position of the 
particle track. 

The helium ions produced by the charged particle and the resulting avalanche create 
UV photons which travel further along the cell. These cause further ionisation in 
the drift cell and new avalanches. Under the high voltage environment of the drift 
cell this process is smooth and propagates at a constant velocity toward both ends 
of the cell. Upon reaching the ends, the ions are collected by the copper cathode 
rings. The arrival time of the plasma at each cathode, which is measured by the 
TDC of each cathode, is used to determine the vertical position of the particle track 
passing through the cell. 

Helium is a low atomic number gas, thus it minimises the energy loss by a charged 
particle passing through the gas. Ethyl alcohol plays the role of the quencher, which 
neutralises helium ions and limits the creation of the UV photons and therefore 
avalanches. A small amount of argon (1%) and water (1500 ppm) were added later 
during the commissioning of the detector in order to increase the plasma propagation 
efficiency and to reduce noise and improve the cell stability [57]. 

4.4 The calorimeter 

The calorimeter in NEMO 3 measures the energy and the time of flight (TOP) of 
particles. It consists of 1940 large scintillator blocks coupled to 3" and 5" photomul- 
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tiplier tubes (PMTs) via light guides. The scintillator blocks are produced from low 
radioactive polystyrene which also provides minimum back scattering. The blocks 
cover the inner and outer cylindrical walls that surround the tracking volume and 
the top and bottom of each sector. In order to minimise energy loss, the scintillator 
blocks are mounted inside the helium-alcohol gas mixture of the tracking volume. 
Figure 4.7 shows a schematic view of one of the scintillation detectors with a 5" 
PMT. The energy resolution (full width half maximum, FWHM), c^Z-E, of the 
NEMO 3 calorimeter is 14% at 1 MeV and 8% at 3 MeV, and the time resolution 
is 250 ps. In addition to 1940 PMTs in the calorimeter, there are six reference 
PMTs situated outside the detector that are only used during laser calibrations 
(Section 4.8). 

4.4.1 The scintillators and the light guides 

The scintillators of the calorimeter are made in seven different shapes designed to 
completely cover the cylindrical geometry of NEMO 3. The thickness of the scintil- 
lators is approximately 10 cm throughout. The width and length are between 11 cm 
and 15 cm. Charged particles lose energy rapidly when passing through polystyrene 
via molecule excitation. As the molecules de-excite, photons are produced at a 
wavelength transparent to polystyrene. Photons incident on the scintillator lose 
energy via Compton scattering and the resulting electrons then lose energy as de- 
scribed previously. The scintillation detectors have a detection efficiency of 50% at 
500 keV. 

In order to transmit scintillation photons to the PMTs, polymethyl methacrylate 
(PMMA) light guides are optically glued to the back face of the scintillator blocks. 
The light guides are 60 mm thick and are glued to the detector walls to provide a 
seal against the helium-alcohol environment of the tracker. 
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Figure 4.7: A schematic view of a scintillator counter with a light guide and 5" PMT 
in the NEMO 3 calorimeter [50]. 
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4.4.2 The photomultiplier tubes 



The photomultiplier tubes (PMTs) convert optical signals from scintillators to mea- 
surable electrical signals. They achieve this via the photoelectric effect and avalanche 
multiplication. There are two sizes of PMTs used in the calorimeter, 3" and 5" . The 
3" PMTs are mounted on scintillators at the top and bottom regions of each sector, 
and the 5" PMTs are mounted on scintillators around the inner and outer cylinder 
walls. The PMTs used have been specially produced to have three times lower ra- 
dioactive contamination than a standard PMT. The gain of the PMTs is adjusted 
to cover a range of energies up to 12 MeV. In order to protect them from the exter- 
nal light the PMTs are contained in black plastic boxes. The output signals from 
PMTs pass through discriminators. These produce logic high signals when the PMT 
pulse passes a minimum threshold. The discriminator signals are used to start the 
time measurement and charge integration which are performed by time to digital 
convertors (TDCs) and analogue to digital convertors (ADCs), respectively. 



4.5 The Trigger and readout 

The trigger ensures that only events of interest are read out. As NEMO 3 is a low 
count rate experiment a simple hardware trigger is used to achieve this. The trigger 
is connected to the calorimeter electronics. When it receives a signal from a PMT 
which has passed a discriminator threshold of 48 mV, it generates the STOP-PMT 
signal after 20 ns. This signal stops the charge integration and time measurements 
performed by an ADC and a TDC, respectively. 

In the next stage, the tracking readout system is programmed to search for activation 
of drift cells in each half-sector (inner or outer). If enough cells are activated (3 out 
of 9 layers), the STOP-A signal is generated and sent to the drift cell acquisition 
boards with a programmable delay set at 6.14 /is after the STOP-PMT. The STOP- 
A signal stops the TDC measurement of the drift cell anodes (TDC^) and the anode 
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drift times are calculated. 

Finally, the STOP-a signal is sent to the drift cell acquisition boards with a fixed 
delay of 710 /is after the STOP-PMT signal. This signal is to stop the TDC^, which 
is a TDC independent of TDCa, and is used for calculation of the drift time for 
delayed hits. This is designed to detect alpha (a) particles from daughter isotopes 
of radon and have a half-life of 164 /is. The value of 710 /xs is more than four times 
this half-life to have high detection efficiency for a-particles. 

4.6 Magnetic coil and shielding 

A solenoid surrounding the entire detector produces a 25 Gauss magnetic field ver- 
tically through the detector. The coil is cylindrical and is 5220 mm in diameter 
and 2713 mm in height. This causes charged particle tracks to bend as they pass 
through the detector allowing the identification of electron and positron tracks. 

Iron plates, wood panelling, borated water and an anti-radon tent form the external 
shielding for the detector. The iron plates are 20 cm thick and use low radioactivity 
iron to stop photons coming from outside the detector. The water and wood com- 
prise a thick shield which slows down fast neutrons and captures them. The wood 
covers the top and bottom of the detector and is 28 cm thick. The borated water 
tanks, which are 35 cm thick, cover the cylindrical external walls of the detector. 

4.7 The anti-radon facility 

After running the detector for approximately one year, it was discovered that the 
radon level inside the detector was too high (with a total activity of around 0.7 Bq) 
and adversely affecting the experiment's sensitivity. This was caused by radon 
diffusion through the glued joints between sectors. In order to reduce the radon 
contamination inside the detector an anti-radon facility was installed in September 
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Figure 4.8: Evolution of the total radon activity in the NEMO 3 before (Phase 1) 
and after (Phase 2) the installation of the anti-radon facility [58]. 

2004. It consists of an airtight tent and a radon trapping facility. 

The detector is fully enclosed in an airtight tent that is made of two layers of 
polyethylene. It allows only radon-free air to pass through from the radon-trapping 
facility into the detector and isolates it from laboratory air. The radon trapping 
system uses activated charcoal that has been treated with oxygen to open up nu- 
merous pores inside the material. When air passes through the charcoal, radon is 
trapped in these pores. The trapping time is greater than the radon decay time. 
Therefore radon decays before reaching the tent. 

As a result of this system, the radon activity inside the detector has been reduced 
by a factor of six. Figure 4.8 shows the level of radon inside NEMO 3 before and 
after the installation of the anti-radon system. 



4.8 Calibration of the calorimeters 

In order to convert the charge output of the PMTs to energy, a calibration procedure 
is performed. In NEMO 3 there are dedicated calibration runs which take place once 



49 



per month and last 24 hours. During this time radioactive sources with well known 
energy spectra are introduced near the double beta decay sources. As NEMO 3 is 
mainly focused on electron energy measurements, the selected radioactive sources 
emit electrons. The ^^''Bi source decays by a conversion process to two electrons [59] 
(see Section 6.1.1) with energies of 482 keV and 976 keV, and is suitable for energy 
calibration up to 1.5 MeV. In order to calibrate the energy up to 3 MeV and higher, 
^°Y, which is a pure beta emitter with transition energy of 2.283 MeV, is used. 

The timing responses of the PMTs are not identical. The time calibration is per- 
formed to make the timing response uniform for all the PMTs in the calorimeter. 
Currently, electrons and photons from ^°^Bi sources are used for the time calibra- 
tion. The particle with the higher energy decays with a delay of 133 /is. This time is 
subtracted from the arrival time of the delayed signal. The time difference between 
the prompt and delayed signals is then used to make the time response uniform. 

Daily studies of the stability of the PMTs (in time and energy measurements) are 
performed using a laser survey system. The purpose of this system is to check the 
absolute energy and time calibration. It also measures the PMT response linearity 
between and 12 MeV. Laser light with a known intensity is sent to the PMTs of 
the calorimeter via optical fibres. Six reference PMTs are connected to scintillators 
with embedded ^°^Bi sources and also to the laser survey system via optical fibres. 
They are used to check the stability of the laser light intensity which is received by 
the PMTs in the calorimeter. The laser energy calibration is not applied to the data 
used in this thesis, as the laser data analysis was not ready for use with the analysis 
presented. The effect of not using this correction has been studied and is included 
as systematic uncertainty (Sections 7.3). 

The laser time correction 

During the analysis of NEMO 3 data in September 2006, it was discovered that 
the timing measurements of the PMTs were not correct. The time measured by 
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the TDC counter of a given PMT was expected to remain constant for different 
laser runs; however, a discontinuity was observed in the distribution of the TDC 
time (Figure 4.9). This effect was seen in the PMTs of all the 20 sectors. To solve 
the timing problem a correction was introduced known as laser time correction 
(LTC) [60]. The measured TDC for each PMT is compared to the reference run as 
follows: 

LTC} = tdc} - td4o4 (4.1) 

where tdcj is the time measured by the PMT i in the laser run j, and 1404 is the 
reference laser run. In this run the maximum number of PMTs were active and they 
were known to be relatively stable and the timing measurements performed by the 
TDCs were known to be correct. The value of LTC is added to the TDC values of 
the PMTs and calculated in each laser run to correct the deviation in timing. 
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Figure 4.9: The mean TDC values for PMTs in outer part of the sector 5 versus 
the laser run number for years 2003, 2004, 2005 and 2006 [61]. The TDC jumps are 
seen from the end of 2003 onwards. 
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Chapter 



5 



Analysis 



technique 



This chapter describes techniques used to prepare data for analysis, including the 
reconstruction of raw data and Monte Carlo (MC) simulation. The ^^°Nd foil ge- 
ometry is defined in Section 5.5. Section 5.6 describes the time of flight (TOF) 
selection that is used in all analysis channels. The statistical methods used to com- 
pare simulated events to data and set limits on several Ou^P processes are described 
in Section 5.7 and Section 5.8, respectively. 

5.1 Reconstruction of particle tracks 

The raw hits in the tracker that are written to data storage tapes are processed by 
reconstruction software. According to the tracker anode and cathode time values, 
the triggered drift cell hits are classified as "in time hits", "delayed hits" or as 
cells fired by neighbouring cells ("re- fired cells" and "noisy cells"). "In time hits" 
are usually due to the electrons or positrons and "delayed hits" are caused by a 
particles. Hits are only reconstructed as tracks if they are "in time hits" or "delayed 
hits" . For "in time hits" the anode time values must not be more than a maximum 
value (currently 1.74 /is found from Geiger calibrations). For "delayed hits" , the time 
measured by the anode is up to 710 fis. The anode drift time and the cathode times 
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provide the transverse and longitudinal position of a hit in a drift cell, respectively. 

Tracks are reconstructed if there are "in time hits" in at least three of the nine drift 
cell layers, with at least two of the hits occurring in neighbouring layers. The pattern 
recognition is carried out using a cellular automation algorithm [62] . This algorithm 
defines segments as the lines which connect two hits in neighbouring layers. In order 
to take into account all inefficiencies, the segments may connect two hits which are 
not in neighbouring layers if there is no hit in between. The angle between the two 
segments with a common hit can not be more than 40 degrees. To connect two cells 
with a segment there are four possibilities, as each cell has two edges, right and left. 
The ambiguity is resolved by finding the longest and smoothest pattern of combined 
segments. 

Once the pattern is found, an iterative fit is performed assuming the track follows a 
helix (due to the magnetic field). The curvature of the track indicates if the particle 
is an electron or positron. The reconstructed track is extrapolated back to the source 
foil and projected onto an associated scintillator surface. The coordinates where the 
track crosses the foil and the scintillator are calculated. Assuming that the tracks 
originate from the foil, the former gives the vertex of the event. For "delayed hits", 
the time which is measured by TDCq, is used to reconstruct the track in the (x, y) 
plane. 

The vertex resolution 

The resolution of the reconstructed vertex of the two electrons in the transverse plane 
depends on their energy. By using two electrons from ^°^Bi source foils (Section 4.8), 
the vertex resolution is found to be 3 mm at 0.5 MeV and 2 mm at 1 MeV in the 
(x, y) plane. The longitudinal vertex resolution depends on both energy and position 
of the electrons in the z direction. For 1 MeV electrons, the resolution is 5 mm if 
the vertex is in the top or bottom of the detector and 4 mm if it is in the centre. 
For 0.5 MeV these values change to 8 mm and 6 mm, respectively [63]. 
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5.2 Particle identification 



Particle identification is possible by combining the information provided by the 
tracker and the calorimeter. NEMO 3 is able to identify electrons, positrons, photons 
and alpha particles. 

Due to the magnetic field charged particle trajectories are curved. In NEMO 3 the 
curvature sign of a trajectory is found with respect to the surface of the source 
foil. Thus, an electron (positron) originating from the foil is identified as a track of 
negative (positive) curvature, which is associated with only one scintillator hit. A 
photon (7) is identified as a scintillator hit which is not associated with a track. As 
the energy losses between the scintillator blocks are not known, events are rejected if 
the photons or electrons are detected by two scintillators. The rate of fake scintillator 
hits is negligible. 

An alpha (a) particle can be identified only if it is accompanied by at least one 
electron track. Alpha particles are mainly produced by ^^^Po decay which is in 
the radon (^^^Rn) decay chain. They are defined by drift cell hits delayed by more 
than 70 /is. As the range of a-particles from ^^^Bi decay in the helium gas of the 
tracking chamber is 25-40 cm [61], the distance between the reconstructed vertex of 
the electron and the delayed drift cell(s) corresponding to an a-particle is required 
to be \^z\ < 30 cm in the vertical direction and |Ar| = a/ (Ax)^ + (A?/)^ < 25 cm 
in the x — y plane. 

5.3 Event simulation in NEMO 3 

Monte Carlo (MC) programs are used as event generators that simulate initial par- 
ticles from isotope decays. The generated events are passed through simulations 
of the detector response so that the output can be compared to data. Resolution 
and detector acceptance effects are taken into consideration in the simulation. MC 
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simulations are necessary for analysing experimental data. The MC events are used 
to measure detector acceptances and to compare the generated physics signal with 
that seen in the detector. Selection criteria can be optimised by studying the signal 
to background ratio for specific processes. 

The NEMO 3 simulation program [64] uses GENBB as the generator. GENBB 
simulates the initial kinematics of the particles for two-neutrino double beta decay, 
different theories of neutrinoless double beta decay and all the possible radioactive 
backgrounds to the double beta decay processes. GENBB also provides the pos- 
sibility to generate the kinematics of Compton scattering of photons from sources 
outside the detector (external photons) and M0ller scattering of external electrons. 

The particles generated with GENBB (photons, electrons, positrons and alphas) 
pass through the various regions of the detector. The description of the geometry 
of the detector and the simulation of the detector response is developed in the 
framework of the GEANT 3 package [65] . The interactions of the particles with the 
source foils, tracking wires, the scintillators and other material in the sectors are 
taken into account in the simulation. 

5.3.1 The reconstruction of the simulated events 

The analysis in this work involves estimating the composition of data samples based 
on MC simulations of the various processes. The simulated events were reconstructed 
in the same way as the experimental data. The functionality of the calorimeter and 
the tracker components in each run period is slightly different. This gives different 
detection efficiency for each run (for example having noisy PMTs disconnected from 
the detector decreases the efficiency). For the reconstructed MC simulation to have 
real detector conditions, the simulated events are assigned to the real data run 
periods. The number of events assigned to each run depends on the duration of 
the run. The detector conditions during the particular run are then applied to the 
simulated events associated to that run. 
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5.4 The data set 



The analysis is performed on data taken between 14* February 2003 and 3P* De- 
cember 2006. This corresponds to runs 1869 to 5468. Each run represents a period 
of time in which data acquisition was performed. The length of the runs is typically 
12 hours. Data are removed from the analysis if any of the following conditions are 
met: 

• There was a major problem with the electronics of the detector during the 
run. 

• The run has been taken less than 24 hours after a general electronics shut 
down, as the PMTs of the calorimeter need time to stabilise. 

• They are taken less than 24 hours after ^^''Bi calibration runs. During such 
calibrations PMTs are known to have a high counting rate and require time 
to stabilise [66]. 

• No laser run has been performed on the day that data were taken. The laser 
runs are vital as they give information about the behaviour of the PMTs, and 
also give the value of LTC (see Section 4.8). 

• No appropriate energy or time calibration of the PMTs is available. 

After removing these data, the effective data taking time is calculated to be 924.7 
days. The data are divided into two sets: Phase 1 and Phase 2. Phase 1 is the data 
taken before installation of the anti-radon facility in September 2004 and Phase 2 is 
the data taken after the installation. The Phase 1 and Phase 2 data are combined 
throughout this thesis unless stated otherwise. 
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5.5 Definition of the ^"Nd source foil boundaries 



To define the foil position in tlie transverse plane, tlie polar angle, 0, is replaced by 
a sector number. The sector number is related to the polar angle by: 

20 

sector number = — ^6. (5.1) 

In the construction design of ^^°Nd, this value is between 5.74 and 5.87. However, by 
studying the distribution of event vertex positions, it is observed that data are shifted 
with respect to the defined MC geometry. Figure 5.1a shows the distribution of the 
vertex position for two electron events coming from the foil (higher sector number). 
The discrepancy between data points and MC simulation can be observed on the 
right side of the ^^°Nd foil. From these plots no conclusion can be drawn for the left 
side of the foil, since two-electron events from ^°°Mo overlap with events from ^^°Nd. 
To check if the same effect is also observable on the left side, the electron-photon 
decay channel was analysed. In this channel the event rate for ^'^'^Mo is lower than 
^^°Nd as this isotope has less contamination than ^^°Nd. Figure 5.1b shows the 
distribution of the vertex position for one electron and one photon events coming 
from the foil (full details of the electron-photon selection are given in Section 6.1.1). 
The discrepancy between simulation and data is observed on both left and right sides 
of the ^^°Nd foil. In order to correct this discrepancy, the ^^°Nd sector number was 
shifted by 0.01 of a sector number to lower values, in the simulation. The effect of 
this shift is shown in Figures 5.2a and 5.2b. The jndj values of the data and MC 
comparison improve significantly. Therefore the ^^''Nd position was redefined to be 
between sector number 5.73 and 5.86. The active height of the ^^°Nd foil is 234 cm, 
thus the boundaries of ^^''Nd in the z coordinate are defined to be l^;! < 117 cm. 
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Figure 5.1: The vertex position of the a) two electron and b) one electron plus 
one photon events from ^^''Nd and part of the ^'^Ca, ^^Zr and ^°°Mo foils. The 
original MC geometry is shown in red and NEMO 3 data are shown as points. The 
statistical uncertainties on the data points are shown with error bars. The shift of 
data is clearly seen in both plots. 
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Figure 5.2: The vertex position of the a) two electron and b) one electron plus one 
photon events from ^^''Nd and part of the ^^Ca, ^^Zr and ^°°Mo foils. The new 
definition of the MC geometry is shown in red and NEMO 3 data are shown as 
points. The statistical uncertainties on the data points are shown with error bars. 
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5.6 Time of flight selection criteria 



It is possible that a particle from a source outside the foil deposits energy in one 
scintillator block, enters the tracking volume, interacts with the source foil, scatters 
off it and is then detected by another scintillator. These types of events, which are 
known as external events, can mimic two electron events originating from the foil 
(internal events). In order to reduce the number of these background events, a time 
of flight analysis is performed [67]. In this analysis, the time differences between the 
two calorimeter hits are calculated for two hypotheses. In the first hypothesis the 
two calorimeter hits and tracks are caused by two electrons from the foil (internal 
hypothesis). The second hypothesis assumes that the hits are from external sources 
(external hypothesis). The time differences calculated for both scenarios are then 
compared with the measured time difference of the calorimeter hits. 



For each hypothesis, Xth defined as: 



2 _ (Atmes - Atffc) 
Xth — 2 T 2 ' \'~'-^) 



where At^es is the measured time difference between the two scintillator hits; Attn is 
the calculated time difference for the internal and external hypothesis; and ames and 
ath are the uncertainty on the measured and calculated time differences, respectively. 
The value of Atth is related to the length of the particle trajectories, li and I2, and 
the energies measured by the calorimeter, Ei and E2. In the case of the internal 
(external) hypothesis for two-electron events, /Sith is defined as: 

AC* = ^ - (5-3) 

Pi P2 

Atr = ^ + (5-4) 

where the relativistic factor, is related to the energy measured by each calorimeter 
block [Ei) and the electron rest mass (mo): 

J^i + mo 
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The theoretical uncertainty, ath, is found from differentiation: 

where tlf^ = j^, as is the uncertainty on the energy measurement and ai is the 
uncertainty on the track length. 

The scenario for internal and external electron-photon events is similar except pho- 
tons do not leave tracks in the tracker. To measure the length of the photon path, 
the distance between the scintillator hit and the intersection point of the electron 
track with the foil (the electron-photon vertex) is found. The length of the photon 
path, and the time calculated for the photon to traverse this length, t^^, are found 
by: 

^7 \/ (^-^ve •^scin)'^ ~l~ (.Vve Uscin)'^ ~l~ (^tie ^scin)'^ i and — , (5.7) 

where x^ei live and z^e are the coordinates of the vertex position of the events from 
the detector centre, and Xscin-, Uscm and Zsdn are the coordinates of the scintillator 
hit by the photon. As it is not known where in the scintillator the photon hits, these 
coordinates correspond to the centre of the face. The uncertainty on tj^ is related 
to the maximum dimensions of the scintillator face and the energy measured by the 
scintillator. For an internal event, Atmes — ^^th* expected to be zero ns with 
some uncertainty due to the time resolution, whereas Atmes — Atl^^ has a value of 
the order -5 ns. The opposite applies for external events. Figures 5.3a and 5.3b show 
Atmes — Atl^^ versus Atmes — Atf^^ for internal two electron events (simulated 2i^/5/5 
decay of ^^"Nd) and external two electron events (simulated ^^^Bi contamination in 
the PMTs), respectively. 

For both hypotheses, the probabilities, P{Xint) -^(xLt)' calculated. The 
P(x^) is defined as: 

P{X') = 1-^ r e^'dx, (5.8) 
V^r Jo 

where x = ^ ^JY^ ' -^^S^^^ ^•'^ shows the P{Xint) ^^'^ ^ixlxt) ^'^^ simulated internal 
two electron events (2i//3/3 decay of ^^°Nd) and simulated crossing electrons events 
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Figure 5.3: Distribution of At^es — AtJ^* versus At^es — At^^* for two tracks from 
a) 21/(3(3 decay of ^^''Nd foil and b) external electrons originating from ^^^Bi contam- 
inants in the PMTs. The simulations are normalised by an arbitrary factor. 
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Figure 5.4: The distribution of the internal a) and external b) hypotheses probability, 
P(xL) and PixLt), for 2u(3(3 events from ^^^Nd foil, c) The P(xL) and d) P(xLt) 
distribution for two electron tracks created by ^^^Bi contamination in the PMTs. 
The simulated events are normalised by arbitrary factors. 

(^^^Bi contaminant of the PMTs). Figure 5.5 shows P{Xint) and P{xlxt) fo^ simu- 
lated internal electron-photon events (^^^Bi contaminant in ^^°Nd foil) and external 
electron-photon events (^^^Bi contaminant of the PMTs). 

By using these MC distributions, data selections are optimised in order to select 
events originating from the foil or from the detector components [68]. The criteria 
for the internal events are 

P{x%nt>OM and P(xL) < 0-01, (5.9) 

and for the external events are 

P(xL)<0.01 and P(xL)>0.04. (5.10) 
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Figure 5.5: The distribution of the internal a) and external b) hypothesis probability, 
PiXint) ^^'^ -P(xLi)' ^7 events from ^^^Bi contaminant inside the ^^°Nd foil, and 
for c) and d) 67 events created by ^^^Bi contamination in the PMTs. The simulated 
events are normalized by arbitrary factors. 

5.7 Fitting Monte Carlo samples to data 

Chapters 6 and 7 are dedicated to the measurement of the radioactive background 
activities and 2z//5/5 half-life by comparing MC simulations to data. In the situation 
where only one MC component is not known, it is sufficient to do a simple fit based 
on the total number of data and simulated events. In this method the simulated 
events are normalised to the same number of events as total data events. 

However, in the situation where more than one MC components should be estimated, 
a binned maximum likelihood fit [69] is used. In this technique, to fit m number 
of MC sources to data, the shapes of the distributions are taken into account. The 
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predicted (normalised) number of MC events in bin i is given by 

f^ = Nnj:^, (5.11) 

j=i J 

where No is the total number of events in the data sample, Nj is the total number 
in the MC sample for source j, Pj is the proportion of source j in data, and Aij is 
the expected number of MC simulated events from source j in bin i (which includes 
the unknown uncertainty in the bin). 

The probabilities of observing a particular number of data events, di, and a particular 
number of simulated events, aji, follow Poisson statistics: 

P{d,-J,)=e-f^^, (5.12) 
dil 

P(a,,;A,,) =e-^-^. (5.13) 
The total logarithm of likelihoods is then defined by: 

n n m 

InL = \nP{di; fi)+lnP{aji;Aji) = ^ c/j In /j - /j + ^ ^ a^i InAjj - Ajj. (5.14) 

i=l 1=1 j=l 

The proportion of each MC source, Pj, is estimated by maximising Equation 5.14 
using computing packages [70, 71]. 



5.8 Limit setting 

The search for new physics in a double beta decay experiment involves dealing with 
small numbers of expected signal events (in the case of OujSP) and much larger 
number of background events than expected signal events (in the case of Oz//5/3 with 
the emission of Majoron(s)). The exclusion of a signal at a particular confidence 
level (CL) can be significantly improved relative to a simple counting experiment 
by using the binned distribution of data, signal and background. In this way each 
bin is treated as an independent search channel. 

A typical new physics search analysis is described by a final variable which is de- 
signed to distinguish between two possible hypotheses: the production and detection 
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of the new physics along with the background (S+B hypothesis), or the presence of 
only the background (B-only hypothesis). The consistency between data and the 
signal and background models in this final variable is used to set a limit on the 
maximum number of signal events consistent with the data. The likelihood ratio 
test statistics [72, 73] is an optimal choice for searches with small statistics [74], 



where L{S + B) and L{B) are the Poisson likelihood for the S+B and B-only hy- 
potheses. 

The limits in this thesis are calculated by using a likelihood-fitter [75] developed at 
the D0 experiment using a log-likelihood ratio (LLR) test statistic. A description 
of this method is given in [76]. The LLR is defined as [72]: 



where i is the the i bin in the distribution used, Si is the number of expected 
signal events in that bin, bi is number of expected background events and Ui is 
the number of data events observed in bin i. The expected LLR distributions for 
both hypotheses are generated by performing a large number of toy Monte Carlo 
pseudo-experiments. In pseudo-experiments, the pseudo-data are the coherent sum 
of signal and background or background processes and thus their value in each bin is 
obtained by drawing a random number from a Poisson distribution where P{x]p) = 
p^Q-p lx\ is the probability for obtaining x events, given an expectation of p events. 
The expectation of a number of events in each bin for each pseudo-experiment is 
varied according to the systematic uncertainties. The systematic uncertainties are 
introduced into signal and background expectations, via Gaussian distributions. 



\ i / 

where is a random number taken from a Gaussian distribution with a mean of 
zero and width of one and (jj is the j^^ fractional uncertainty on the rate of the . 



Q = L{S + B)/L{B), 



(5.15) 




(5.16) 




(5.17) 
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The same random Gaussian number is used for correlated uncertainties. The LLR 
distributions for S+B and B-only hypotheses, buih up from the pseudo-experiments 
can be compared with the observed value of LLR in the data (LLRobs or Xd)- The 
LLRobs are found by substituting Ui in Equation 5.16 with the number of observed 
data events (rfj). The confidence level in S+B hypothesis is given by 

CLs+B = PsMx>Xd)= -4^dx, (5.18) 

Jxd dx 

where Ps+b{x > Xd) is the probability for the S+B hypothesis to produce an out- 
come which is less signal-like than that observed in the data and is defined by 
integrating the expected LLR distribution for S+B hypothesis. Similarly, the cal- 
culation of the confidence limit for the background only hypothesis is given by the 
probability for the B hypothesis to produce an outcome which is less signal-like than 
that observed in the data and is found by integrating the expected LLR distribution 
for B-only hypothesis, 

CLb = Pb{x >Xd)= -T^dx. (5.19) 

Jxd dx 

Downward fluctuations of the background lead to inconsistencies between expected 
background and observed data and create inaccurate exclusion limits using CLs+b [74]. 
To deal with this problem, CLs+b is divided by CLb and the value of the signal 
confidence level (CLg) is defined as [72]: 

CLs = (5.20) 

ULb 

Figures 5.6a and 5.6b show two examples of LLR distributions. The data originating 
from ^^°Nd foil are used to create these distributions. In both figures the background 
model consists of 2z//5/5 plus all other radioactive background (See Chapter 7). In 
Figures 5.6a, the signal is Oh'f3(3 and in Figure 5.6b, it is 0^(3(3 with the emission of 
a Majoron particle. 

In order to calculate the observed limit on the signal rate, this rate is given a new 
value and the expected LLR distributions are re-generated from another series of 
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pseudo-experiments and the observed LLR is recalculated. This leads to a new 
value of CLg. This process is repeated through a range of signal rates until CL^ 
reaches the desired confidence level, which in this thesis, by convention, is CL^ = 0.1 
(1 - CL, = 0.9) or 90% confidence level. 

5.8.1 Definition of observed and expected limit 

There are two distinct limits, the median expected and observed limit. The values 
of observed LLR (LLRobs) and median expected LLR (LLRmcd) are found by setting 
Hi = di and rii = hi in Equation 5.16, respectively. The value of LLRobs relative to 
LLRmed indicates whether data appears to be more background-like or not. The ex- 
pected limit on the signal rate gives the experiment sensitivity without any bias from 
actual observation. Thus, different experiments searching for a particular process 
can be compared using the expected limit results. 

5.8.2 The profile likelihood ratio 

Systematic uncertainties are often estimated rather than measured and frequently 
correspond to an upper bound on the true value. These systematic uncertainties 
widen the expected LLR distributions and decrease the analysis sensitivity for set- 
ting limits. This effect of systematics can be reduced by comparing the predicted 
shape of the final variable to that observed in data, and finding the best fit for 
systematic uncertainties to data observation [77]: 



where p{H)[ is the systematically varied prediction in bin i for the two hypotheses: 

k 

p{H)\=p{H),\[{l + S{H)u). (5.22) 

k=l 

In the case of the S+B hypothesis Pi = Si + bi, and in the case of the B-only 
hypothesis Pi = bi. The fitted value for systematics k is indicated as S{H)k. The 




k 



(5.21) 
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Figure 5.6: LLR distributions for the expectation in the B-only hypothesis (blue 
hne) and in the S+B hypothesis (red hne) and the observation in the data (green 
hne) for a) OujSf] and b) 0i'j3(3x search. The regions integrated to obtain 1 — CLb 
and CLs+B are also shown. 



70 



optimal choice of the test statistic in this case is the profile likelihood ratio [75, 78]. 
This test statistic is the ratio of the minimised for each of the two hypothesis. The 
LLR is defined as: 

LLR=-ln fr"[f^\ , (5.23) 

To find the LLR distributions for the S+B and B-only hypotheses, rij in Equa- 
tion 5.21 is replaced with the MC pseudo-data, and to find LLRobs and LLRmcd, 
rij is replaced with the number of observed data and expected background events. 
The p{H)^ term in Equation 5.21 results in narrower LLR distributions than the 
ones resulting from Equation 5.16 and therefore increases the analysis sensitivity for 
setting limits. The Equations 5.18 to 5.20 are used to set limits on the number of 
observed and expected events at 1 — CLg = 0.9 by varying the signal rate as before. 



5.8.3 Floating the background normalisation 

The function used in the fit model in Equation 5.21 includes a term {S{H)l) that 
adds a Gaussian constraint on each systematic uncertainty . If the normalisation on 
the background is not fully known, it can be "floated" by removing the Gaussian 
constraint. This allows a free normalization of the background in the fit. 
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Chapter 6 

Estimation of the radioactive 
background 



This chapter describes the measurement of the radioactive backgrounds to double 
beta decay of ^^°Nd in NEMO 3. The dominant decay channels of each background 
are discussed. The backgrounds to OuPP are divided into 2z//5/5 background and 
radioactive backgrounds that can mimic double beta decay. The former is discussed 
in Chapter 7. The latter, which are also backgrounds to 2i'j3f3, are divided into two 
categories: internal and external. The internal backgrounds are due to radioactive 
contamination inside the NEMO 3 source foils, whereas the external background is 
defined as events created by the radioactivity located in the detector's components 
and the material surrounding the detector. Most contamination comes from the 
two main natural decay chains of uranium (^^'^U) and thorium (^^^Th) . These decay 
chains are shown in Figure 6.1. The detector components and the source foils are also 
contaminated with potassium-40 (^''K), which is a naturally occurring radionuclide. 



6.1 The internal background of the Nd foil 

The internal contaminants can mimic double beta decay of ^^°Nd in four ways: 
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Figure 6.1: The three main decay chains of natural radioactivity (uranium-238, 
thorium-232 and uranium-235) [50]. The two former are the main source of back- 
ground to double beta decay in NEMO 3. The isotopes in grey produce the main 
background to neutrinoless double beta decay. 

• A /5 decay accompanied by an electron from an electron conversion process. 
The latter happens when a photon from an excited daughter nucleus interacts 
with an electron in one of the inner electron shells, causing the electron to be 
ejected from the atom. 

• M0ller scattering of a /3-decay electron in the foil which leads to emission of 
two electrons. 

• A /3-decay accompanied by a photon. In this case a second electron is produced 
via a Compton scattering. 

• Two-electron conversion processes due to the de-excitation of the daughter 
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Figure 6.2: Four ways that the internal contaminants inside Nd source foil can 
mimic double beta decay. 

nuclei of ^°^Bi, ^^^Eu and ^^^Eu contaminants in the ^^"Nd foil. 

Figure 6.2 shows the schematic of each of these four processes that can give two 
electrons with sufficient energy from a single vertex. 

The internal contaminants vary in the different source foils of NEMO 3. The level 
of contamination was obtained (after the isotope purification) from measurements 
using a low-background high purity germanium (HPGe) detector [50]. This HPGe 
detector was able to detect photons from de-excitation of the daughter isotopes. The 
energy and the intensity of the photons indicated the contaminant and its abundance 
in the foil. 

The HPGe radioactivity measurements for the Nd203 source foil are shown in Ta- 
ble 6.1. The half-life (T1/2) and the nuclear transition energy (Q) are also included 
for each element. All of these contaminants are f3 emitters, except ^^''Bi. The ^^^Bi 
and 208x1 with high Q value are the main radioactive background to OujSlS of ^^°Nd. 
The other elements given in Table 6.1 can mimic 2uj3l3 events, except ^^^Th and 
2i2p]3 whose Q value is not high enough to mimic 2v(3(3. However, ^^^Th decays to 
234m which has a half-life of 1.17 min and undergoes a (3 decay with Q value of 
2.29 MeV and is thus a background to 2u[3(5. The same statement is true for ^^^Bi 
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Table 6.1: Radioactivity measurements of Nd203 source foils performed by the 
HPGe detector, after purification [50]. 

which is a daughter isotope of ^^^Pb and has Q value of 2.25 MeV and half-life of 
60.5 min. 

The HPGe measurements show that the ^^°Nd foil is contaminated with ^°^Bi and 
^^^Eu atoms. Both of these isotopes are artificial as they are not found in nature. 
Thus it is likely that the equipment used for ^^°Nd enrichment was polluted with 
small amounts of this isotope. Europium-152 is created from natural europium and 
is usually produced together with ^^^Eu (natural europium is a mixture of ^^^Eu and 
^^■^Eu) [79, 80]. This means, in addition to ^^^Eu atoms, there should be some ^^^Eu 
contamination, which the HPGe detector was not able to measure. Europium- 
154 has very similar decay schemes and half-life as ^^^Eu. {Q = 1.97 MeV and 
Ti/2=8.6 y). Both isotopes can undergo either P decay (branching ratio of 27%) or 
an electron capture (73%). 

The HPGe detector measurements are not definitive and are performed for a small 
sample of ^^°Nd and before installing the foil inside the NEMO 3 detector. In 
order to have a correct estimation of each internal radioactive background, precise 
measurements of the activities with NEMO 3 data need to be performed. 

Background contaminants often decay in multiple channels (and not only the two- 
electron channel). It is vital to measure the background activities in channels which 
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are independent of the 2z//5/5 signal. The activity of each background is determined 
using the decay channel with the highest branching ratio. The activity (A) of each 
contaminant (in Becquerel) is defined by the number of decays per second, 

^ ^ {Ndata - Nbgr) 

where N^ata is the number of data events found for the decay channel, Nbgr is the 
number of background events to a specific contaminant, t is the total data taking 
period in seconds and e is the efficiency of the event selection criteria, which is 
found by the ratio of simulated events that pass the selections to the total number 
of generated events. 



6.1.1 Electron-photon channel 

The isotopes ^°^Bi, ^^^Eu and ^^^Eu decay predominantly to an electron and a pho- 
ton. ^°^Bi can undergo electron capture and decay to the excited state of ^°^Pb. 
In this case the photon is produced by the de-excitation of the excited ^°^Pb and 
the electron is produced by the internal conversion process [59]. Figure 6.3 shows a 
decay scheme of ^°''Bi, in which only the most intense transitions are shown. This 
decay scheme shows that the electron and photon energy peaks are at 0.57 MeV and 
1.06 MeV. In the NEMO 3 experiment, due to energy loss in the foil, the energies 
of these peaks are slightly shifted. 

^^^Eu can either undergo a (3 decay (27%) or electron capture (73%) to produce 
several excited states of ^^^Gd or ^^^Sm. These daughter isotopes de-excite to the 
ground states resulting in the production of photons of various energies [81]. Simi- 
larly, ^^^Eu undergoes P decay to produce ^^''Gd [79] . 

208rpj (decays predominantly into one electron and two photons; however, it is possible 
that only one of the photons is detected and thus the activity of this isotope can be 
checked using the electron-photon channel. Further details of this isotope's decay 
scheme are given in Section 6.1.2. 
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0.57 MeV 
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Figure 6.3: A decay scheme of ^''^Bi. 

To select electron-photon (67) events originating from the foil the following event 
selection criteria are applied: 

• Only one electron track associated to a scintillator hit with deposited energy 
greater than 0.2 MeV. 

• The vertex of the 67 event, which is found by extrapolating the electron track 
to the foil, is reconstructed inside the ^^'^Nd foil (defined in Section 5.5). 

• The track has at least one hit in one of the first two layers of drift cells. 

• Only one scintillator block is found which is not associated with a track, with 
energy greater than 0.2 MeV. 

• The internal TOF hypothesis is greater than 4%, and the external TOF hy- 
pothesis is less than 1%. 

• The energy sum of all other scintillator deposits is less than 0.15 MeV. This 
requirement rejects events with more than one photon. 

• The TOF calculation is less precise for particles with short trajectories. Thus 
the events with photon and electron trajectory lengths less than 50 cm are 
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Figure 6.4: a) The track length of the electron and b) trajectory length of the photon 
in 67 events. The statistical uncertainties on the data points are shown with error 
bars. 

removed. The photons and electrons with trajectory length less than 50 cm 
are detected by the scintillators situated on top and bottom of the detector. 
Figures 6.4a and 6.4b show the electron track length and photon trajectory 
length before applying this requirement, respectively. Photons can have longer 
trajectory lengths than electrons because they can cross boundaries of sector 5 
and be detected by the scintillators in the neighbouring sectors. 

Figure 6.5 shows the reconstructed z and ^ components of the vertex in the ^^'^Nd 
foil for the 67 events that pass the event selection criteria. Several hot-spots are 
seen in this figure. The two main hot-spots are in the regions: 



1.82 < < 1.827 rad and 34 < ^ < 42 cm, 
1.836 <d)< 1.842 rad and 6 < 2 < 12 cm. 



(6.2) 



Figures 6.6a and 6.6b show the electron and photon energies inside the hot-spot 
region. The normalisations of the ^°''Bi energy distributions are found by fitting 
simulated ^°^Bi events to data after background subtraction. The ^^''Bi activity is 
found to be ten times higher in these hot-spots and therefore this leads to the con- 
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Figure 6.5: Distribution of event vertices in the ^^'^Nd source foil for 67 events. 
Hot-spots are seen in the plot. 

elusion that these hot-spot areas are contaminated with ^°''Bi and they are therefore 
removed from the subsequent analysis. 

Figure 6.7 shows the cosine of the angle between the electron and the photon. A 
single electron from a (3 decay can be scattered by a scintillator block and create a 
photon which is detected by the neighbouring scintillator. Thus events are rejected 
if the cosine of the angle between the photon and the electron is more than 0.9. An 
67 event originating from the ^^°Nd foil is displayed in Figure 6.8. The intermediate 
circle illustrates the source foil and the blue hits illustrate the electron track. The 
scintillator hits are shown in red. 

^°^T1 activity measurement 

The activity of ^o^Tl is found in the energy region in which it dominates over all other 
background. Figure 6.9 shows the photon energy for simulated ^os^jj (normalised 
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a) 




Figure 6.6: The energy of a) the electron and b) photon for 67 events originating 
from the hot-spot region. The components of total background to ^''^Bi (shown in 
grey) are normalised to activities found from 67 and other decay channels. The 
statistical uncertainties on the data points are shown with error bars. 

to an arbitrary number) and the sum of all other internal contaminants. This plot 
demonstrates that in order to measure the activity of ^os^j ^^le electron and 
photon channel, the energy of the photon is required to be more than 2.0 MeV 
to increase signal over background. As the normalisation factors of other internal 
backgrounds at this stage of the analysis are not known this cut is applied in order 
to reduce internal background events from other contaminants. 

After applying the photon energy requirement, the only background to ^os^l which 
remains in the data sample is the radon from the tracker (described in Section 6.2.1). 
The average activity of this background, which was measured by studying the 
electron-alpha channel [58], is 0.45 ± 0.07 mBq. Figures 6.10a and 6.10b show 
the photon and the electron energies for data and the simulated events with photon 
energy more than 2 MeV. In these figures, the distributions for simulated ^^^Tl and 
for data minus radon background are normalised to the same number of events. The 
event selection efficiency for ^^^T\ MC is (0.390 ± 0.004 (stat))%. The number of 
data events that pass the event selection criteria is 193 of which 45.1 ± 9.4 events 
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Figure 6.7: The cosine of the angle between the electron and the photon in the 
67 channel. The peak in cosine of the angle more than 0.9 is mainly due to the 
scattered single-electron events which are shown in dashed green. The statistical 
uncertainties on the data points are shown with error bars. 



NEM03 - RUN: 5136 EVENT: 56206 (2006-90-03 @ 17:30:12) 




Figure 6.8: A top and side view of an 67 event originating from ^^°Nd foil in Sector 5. 
This event is from data taken in August 2006. 
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Figure 6.9: The energy of the photon originating from ^os'jj jg gj^own in cyan and 
the sum of all other background photon energies are shown in dark blue. 

are estimated to be radon background. This gives a ^os^l activity of 

°^T1) = 0.47 ± 0.05 (stat) mBq, (6.3) 

Since the ^^°Nd source foil has a mass 50.7 gr, this value corresponds to 

A' = A/m = 9.27 ± 1.02 (stat) mBq/kg. (6.4) 

The isotopes ^°^T1, ^^^Bi and ^^^Ac are all in the decay chain of ^^^Th. ^^^Bi can 
emit an a particle and decay to ^os^j with a branching ratio of 36% [82]. Thus 
the activities of ^^^Ac and ^^^Bi are 36% of the ^os^pj activity. The systematic 
uncertainty on this measurement and the final results of ^osxi activity will be given 
in Section 6.1.2. 

^°^Bi activity measurement 

To find the ^''^Bi activity, the following energy selections are applied 

0.7 <Ee< 1.1 MeV, 
0.25 < E^< 0.6 MeV. 
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Figure 6.10: The distributions of the a) electron energy and b) photon energy for 
the events with photon energy more than 2 MeV. The error bars show the statistical 
uncertainties on the data points. 

In these energy regions, ^°^Bi dominates over all other background. Figures 6.11a 
and 6.11b show the electron and photon energies in these energy regions. The simu- 
lated ^°^Bi is normalised to data minus the total number of other background events. 
These backgrounds to ^°^Bi are normalised to the activities found by studying dif- 
ferent energy regions (^^'^Tl) or other decay channels (Sections 6.1.3 and 6.2). The 
number of candidate events is 3734. The number of estimated background events 
is 1217 ± 31 (stat). The efficiency is equal to (0.271 ± 0.001 (stat))%. Thus the 
activity of ^°'^Bi is estimated to be 



A(20^Bi) = 11.7 ±0.3 (stat) mBq, 
A' = 230.8 ± 6.2 (stat) mBq/kg. 



(6.7) 
(6.8) 



^^^Eu and ^^"^Eu activity measurements 

The activities of ^^^Eu and ^^^Eu are found by fitting their simulated photon energy 
spectrum to data minus all other background events simultaneously. A binned 
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Figure 6.11: Distribution of the energy of the a) electron and b) photon energy for 
energies 0.7 < -Eg < 1.1 MeV and 0.7 < E^. < 1.1 MeV. The error bars show the 
statistical uncertainties on the data points. 

maximum likelihood fit as described in Section 5.7 is applied. The fit is applied to 
the full range of the photon energy distribution. The total number of 67 data events 
is 14230. The ^°^Bi and ^os^jj activities are normalised to the activities given in 
Equations 6.3 and 6.8. The total number of background events is 12752 ±69 (stat). 
By applying the fit, the number of 67 events from ^^^Eu and ^^^Eu are found to 
be 834 ± 37 (stat) and 643 ± 35 (stat) and the efficiencies for selecting 67 events 
from these isotopes' decays are (0.380 ±0.005 (stat))% and (0.630 ±0.005 (stat))%, 
respectively. The errors on the number of ^^^Eu and ^^^Eu events are determined 
by the fit. The activities of ^^^Eu and ^^^Eu are: 



AC^^Eu) = 2.71 ±0.12 (stat) mBq 
A' = 53.5 ± 2.4 (stat) mBq/kg, 

AC^^Eu) = 1.26 ±0.07 (stat) mBq 
A' = 24.8 ± 1.4 (stat) mBq/kg. 



(6.9) 
(6.10) 

(6.11) 



Figure 6.12 shows the energy of the electron, photon and the sum of the energies 
of these two particles. The simulated ^''^Tl, ^°^Bi, ^^^Eu and ^^^Eu are normalised 
to the activities found in this section. Other background's energy distributions are 
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E,+E,[MeV] 

Figure 6.12: Distribution of the energy of the a) photon, b) electron and c) electron 
plus photon for the events that pass the event selection criteria. The statistical 
uncertainties on the data points are shown with error bars. 

normalised to activities found from other decay channels. These plots demonstrate 
that background MC fits data well in the internal 67 channel. 

6.1.2 Electron plus two photon channel 

The isotope ^os^jj ^^^^ undergo a /? decay to an excited state of ^°'^Pb which de- 
excites to the ground state via two-photon emission. The first photon has an energy 
of 2.6 MeV and is accompanied by a second photon with an energy of 0.58 MeV 
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or 0.86 MeV. Figure 6.13 shows the decay scheme of ^°^TL The main internal 
background to ^os^l in this channel is ^^^Bi which can undergo a (3 decay and form 
an excited state of ^°^Po which de-excites either by emission of two photons with 
energies of 0.61 MeV and 1.12 MeV or a photon with energy 1.76 MeV (a decay 
scheme is shown in Figure 6.14). 

To select one electron plus two photons originating from the foil, the following event 
selection criteria are apphed: 

• Only one electron track associated to a scintillator hit with deposited energy 
greater than 0.2 MeV. 

• The vertex of the 677 event , which is found by extrapolating the electron track 
to the foil, is reconstructed inside the ^^°Nd foil (defined in Section 5.5). 

• The track has a hit in one of the first two layers of the tracker. 




' — ; 2.61 MeV 

171 

' 

208pj^ 

Figure 6.13: A decay scheme of ^osxi showing the relevant energy levels which can 
create 677 events. 
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Figure 6.14: A simplified decay scheme of ^^'^Bi sfiowing the relevant energy levels 
which can create 677 and 67 events. 

• Two scintillator hits are found with the deposited energies greater than 1.8 MeV 
for the higher energy photon and greater than 0.35 MeV for the lower energy 
photon. These energy requirements which are based on the decay scheme of 
^°^T1 are effective to reject internal ^^^Bi background. Figure 6.15 compares 
the distribution of the higher energy photons for the simulated ^^^Bi and ^°^T1. 

• The internal TOF hypothesis for each 67 pair is more than 4%, and the external 
TOF hypothesis is less than 1%. 

• The energy sum of all other photon clusters is less than 0.15 MeV. 

• The length of the electron and photon trajectories is more than 50 cm. 

Figure 6.16 shows a typical internal 677 event display. Figure 6.17 shows the energy 
distributions of the photons and the electron for 677 events that pass the event 
selection criteria. The energy sum of all the decay particles is also shown. Figure 6.18 
shows the cosine of the angle between particles. The simulated ^os^jj jg normalised 
to the same number of events as data minus ^^"^Bi background. These plots show 
that the background MC describes data well in this channel. The expected number 
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Figure 6.15: Energy distribution of the higher energy photon for simulated ^OBrp^ 
and ^^^Bi events. 

of ^^^Bi events is 0.22 ±0.02 (stat). This value is found by normalising the simulated 
^^^Bi events to the activity of 0.169 mBq, which was measured by studying ea events 
originating from the foil [58]. The number of data events is found to be 251. The 
^°^T1 event selection efficiency for this channel is (0.55 ± 0.01(stat))%. The activity 
of 2°^T1 is thus found to be: 

A = 0.57 ± 0.04 (stat) mBq (6.12) 
A' = 11.04 ± 0.79 (stat) mBq/kg. (6.13) 



Measurement of the background uncertainties by comparing the e'j and 
677 results 

To estimate the systematic uncertainty on the measurement of the ^os^jj activity, 
results from the ej and 677 channels are compared. The mean value of the two 
activities given in Equations 6.3 and 6.12 is obtained. The systematic uncertainty 
is found by calculating the deviation of each activity from the mean of the two 
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NEM03 - RUN: 4778 EVENT: 139828 (2006-41-06 @ 02:27:48) 




Figure 6.16: An 677 event display. A top and side view of the event are shown. 
This event is from data taken in May 2006. 

extracted activities. The final ^^^Tl activity is thus measured to be: 



^(208t1) = 0.51 ± 0.03 (stat) ± 0.04 (syst) mBq, (6.14) 

A = 10.06 ± 0.59 (stat) ± 0.79 (syst) mBq/kg. (6.15) 

By normalising ^^^Tl background to this value and refitting, the ^°^Bi, ^^^Eu and 
^^''Eu activities are remeasured to be: 

A(207Bi) = 11.3 ±0.3 (stat), (6.16) 

AC^^Eu) = 2.6 ±0.1 (stat), (6.17) 

AC^^Eu) = 1.24 ±0.07 (stat). (6.18) 



As with ^'^^Tl, the uncertainties on the ^°''Bi, ^^^Eu and ^^^Eu activities are defined 
by the difference from the mean values. Thus the final activity results for these 
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Figure 6.17: The energy distribution of the a) higher energy photon, b) lower energy 
photon, c) electron and d) sum of all particles for 677 events. The error bars show 
the statistical uncertainties on the data points. 
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Figure 6.18: The cosine of the angle between a) the higher energy photon and the 
electron and b) the two photons for 677 events. The error bars show the statistical 
uncertainties on the data points. 
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isotopes are 

A(20^Bi) = 11.5 ± 0.2 (stat) ± 0.2 (syst) inBq, 

A' = 226.8 ± 3.9 (stat) ± 3.9 (syst) mBq/kg, (6.19) 

AC^^Eu) = 2.6 ±0.1 (stat) ±0.1 (syst) mBq 

A' = 51.3 ± 1.4 (stat) ± 2.0 (syst) mBq/kg, (6.20) 

AC^^Eu) = 1.25 ± 0.05 (stat) ± 0.02(syst) mBq 

A' = 24.6 ± 1.0 (stat) ± 0.4 (syst) mBq/kg. (6.21) 

The reconstruction of an additional photon in the 677 channel gives the system- 
atic uncertainty on the photon detection efficiency. Other systematic uncertainties, 
including the uncertainties on the electron detection efficiency (section 7.3), are 
negligible for the purpose of these results. 



6.1.3 Single-electron decay channel 

The isotopes ^^^™Pa and ""^K predominantly undergo (3 decay. Thus the activities 
of these internal contaminants are measured by studying the single-electron (le) 
channel. The following requirements are applied to select le events: 

• Only one electron is found in the event with energy greater than 0.5 MeV. 
This energy requirement eliminates a large amount of low energy events which 
are not relevant for this analysis. 

• The electron track must originate from the ^^°Nd foil. 

• The length of the track is greater than 50 cm. 

• The track has at least a hit in one of the first two layers of the tracker. 

Figure 6.19 shows the reconstructed (j) component versus the z component for single- 
electron events. A hot-spot is observed in the region: 

1.815 <4>< 1.827 and 104 < 2 < 110 cm, (6.22) 
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which is highly contaminated with ^^^'"Pa (Figure 6.20) and has been subsequently 
removed from the analysis. Figure 6.21 shows a (3 decay event display. ^^°Bi from 
the radon decay chain can settle on the tracker wires and undergo a (5 decay which 
can mimic le events originating from the foil. The average activity of ^^''Bi on the 
surface of the wires is 12.61 mBq [58] (more information about this background is 
given in Section 6.2.1). Other internal backgrounds are normalised by the activities 
found in Section 6.1.1. The activities of and ^^^'"Pa are found by fitting the 
data simultaneously with the MC contributions of these two contaminants. 

Figure 6.22 shows the energy distribution of the single-electron events that pass the 
event selections. This figure demonstrates that the background MC describes data 
well in this channel. The activities given from Equation 6.15 to Equation 6.18 are 
used to normalise internal backgrounds to the le channel. As before, the radon 
activity is normalised to average value of 0.45 ± 0.07 mBq (the external background 
components and their activities are given in Section 6.2). Table 6.2 shows the 
number of data events {Nexp)-, number of background events remaining in the data 
found by MC simulations (Nhgr), number of events originating from and ^'^^™'Pa 
found from the fit (A'^4ok and A'^234mpa) and the event selection efficiency. From this 



Contaminant 


^exp 


Nbgr 




Efficiency 


40K 


392123 


179922 ± 698 


161675 ± 569 


18.7% 


234m 


392123 


179922 ± 698 


50525 ±431 


26.1% 



Table 6.2: The internal contaminants of the ^^''Nd foil that mainly decay to an 
electron, the number of expected events, the total number of background events and 
the efficiency of the event selection criteria. The errors on the number of and 
234mpg gygj^^g gj-Q determined by the fit. 
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Figure 6.19: The reconstructed component versus z component for single-electron 
events. 
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Figure 6.20: The distribution of the electron energy in the region 1.815 < (p < 
1.827 and 104 < 2 < 110 cm for Phase 2 data. The statistical uncertainties on the 
data points are shown with error bars. This plot shows that the hot-spot region is 
highly polluted with ^^''"^Pa. 
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NEM03 - RUN: 5008 EVENT: 9683 (2006-70-09 @ 14:23:15) 
5 




Figure 6.21: A (3 decay event display. A top and side view of the event are shown. 
This event is from data taken in July 2006. 




Figure 6.22: Distribution of the energy of the electron for the single-electron events 
that pass the event selection criteria shown in a) linear scale and b) log scale. A 
simultaneous data and MC fit is performed for '^^K and ^^^'"Pa. The error bars show 
the statistical uncertainties on the data points. 
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table, the activities of and ^^^"^Pa are found to be: 

A(^°K) = 10.8 ± 0.1 (stat) ± 0.6 (syst) niBq (6.23) 

A' = 213 ± 2 (stat) ± 12 (syst) niBq/kg, (6.24) 

j^(234mp^^ = 2.42 ± 0.04 (stat) ± 0.14 (syst) niBq (6.25) 

A' = 47.7 ± 0.8 (stat) ± 2.9 (syst) niBq/kg. (6.26) 

The systematic uncertainties on and ^^^"^Pa activities are found by varying the 
energy selections between 0.5 MeV and 1.5 MeV and refitting ^^^Pa and MC 
distributions to data minus other backgrounds after each energy selection. The 
systematic uncertainties are found to be 5% for and 6% for ^'^^™Pa. 



6.2 The external background 

The external background in the NEMO 3 experiment is caused by electrons and 
photons generated outside the source foil. These particles can interact with the foil 
and mimic two electron tracks. An external photon's interaction with the foil can 
produce an electron-positron pair, which can be misidentified as two electrons. A 
photon can also undergo Compton scattering and produce an electron and another 
photon. The second electron can be produced by either M0ller scattering of the 
electron or Compton scattering of the second photon. Figure 6.23 shows these three 
main ways that an external photon can mimic two -electron events originating from 
the foil. 

The second source of external background are crossing electrons which are produced 
by Compton scattering of the external photons within the detector components. 
These electrons can be scattered by the foil and thus detected by two different 
scintillators. As described in Section 5.6, these events can be reduced by the TOF 
requirements. 

The external background is divided into two main categories based on the origin of 
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Figure 6.23: Three ways that the external background can mimic double beta decay: 
the electron-positron pair production due to the photon interaction with the foil (the 
diagram on the left); double Compton scattering (the diagram in the centre); and 
Compton followed by M0ller scattering (the diagram on the right). 



the events: the radon background inside the tracker and the background originating 
from outside the tracker (i.e. the calorimetry, the shielding and the air surrounding 
the detector). 



6.2.1 The radon background 

Radon (^^^Rn) is out-gassed into the air from the surrounding rocks and can enter 
the detector either through leaks between the sectors, or through gas piping joins. 
Bismuth-214 (^^^Bi) is one of the descendants of ^^^Rn which can settle on the surface 
of the tracker's wires. The ^^^Bi undergoes a /? decay to ^^^Po. The latter has a half- 
life of 161 /is and decays to ^^^Pb via emission of an a particle. Thus the activity 
of ^'^^Bi on the wires was measured by studying electron-alpha (ea) events [58] . 

Another source of background which is part of the radon decay chain originating 
from the surface of the wires is ^^'^Bi. This isotope is a /? emitter with Q value 
of 1.1 MeV. Therefore few ^^'^Bi events can mimic 2z//5/5 decay. However, since it 
contributes to the internal single-electron decay channel, it is important to estimate 
its activity. The activity of ^'^^Bi was measured by studying single-electron events 
originating from the surface of the wires [58]. Table 6.3 gives the activity of the 
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Contaminant 




Activity 


(niBq) 






inner foil side 


outer foil side 




Phase 1 


Phase 2 


Phase 1 


Phase 2 


214Bi 


768 ± 10 


177 ±4 


635 ±8 


145 ±3 


210Bi 


15.5 ± 1.2 


14.7 ± 1.0 


13.4 ±0.9 


12 ± 1 



Table 6.3: The activities of ^^^Bi and ^^°Bi originating from the closest layers to the 
foil in sector 5. Due to the installation of the anti-radon facility, the level of the 
^^^Bi on the wires is reduced dramatically for the Phase 2 period [58]. The errors 
given are due to the statistical uncertainties. 

^^^Bi and ^^°Bi originating from the first layer of the drift cells situated in sector 5. 

The other main source of external background is the natural radioactivity of the 
detector components. The PMTs represent an important source of background. 
The HPGe measurements [50] show that the glass of the PMTs are contaminated 
with 2i4Bi, 40^^ 208rjj 228^^^ ^ g^^^y ^j^g external 67 and crossing electron 

events [58] reveals that the iron shield and iron in the sector walls are also contam- 
inated with these isotopes. It is also important to include ^''Co and ^^^'"Pa to the 
external background model [58] in order to reproduce the crossing electron data. 
The magnetic shield of the PMTs and the wall of the sectors are contaminated 
with these isotopes. In addition, before the installation of the anti-radon facility, a 
fraction of external ^^^Bi and ^os^j background originated from the air surrounding 
the detector. Table 6.4 summarises the external background sources and the values 
measured in [58]. 



6.3 Validation of the external background model 

The results presented in Section 6.2 have been obtained by studying foil crossing 
events throughout the detector. In order to validate the external background model 
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External background 




Activity (Bq) 










208 




60Co 


Glass of the PMTs 


324 


72 


27 


1078 


— 


Plastic scintillators 


— 


— 


— 


21.5 


— 


Magnetic shield of the PMTs 


— 


— 


— 


— 


14.6 


Sector iron walls 


9.1 


8.5 


3.1 


100 


6.1 


Internal tower 










18.4 


Iron shield 


7359 


1345 


484 






Air between the detector 


PI: 


566.5 




11.5 






and the iron shield 


P2: 















Table 6.4: Components of the external background model [58]. 

for the ^^''Nd foil, only particles interacting with this foil are selected. Two decay 
channels are used to check the external background model: the crossing electron 
channel and the external 67 channel. In both channels the simulated external back- 
grounds are normalised to the activities given in Table 6.4. The distribution of the 
sum of external background is then compared with data. 

The one crossing electron events create two tracks in the tracking chamber and two 
scintillator hits in the calorimeter. The external TOF hypothesis for these events 
must be more than 4% and the internal TOF hypothesis must be less than 1%. 
Figures 6.24 shows the energy of the crossing electron after scattering by the foil 
and detection by the second scintillator, the sum of the energies deposited in the two 
scintillators and the cosine of the angle between the two electron tracks. Phase 1 
and Phase 2 data are shown separately. The total number of background events 
estimated by MC simulations is 8177 ± 103 (stat) for Phase 1 and 8197 ± 101 (stat) 
for Phase 2. The number of data events are 8050 and 8012. The difference between 
the number of expected background events and data events gives 2.2% uncertainty 
on external background in this channel. 
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Figure 6.24: a, b) energy of the crossing electron after scattering by the foil; c, 
d) sum of the energies deposited in two scintillators and e, f) cosine of the angle 
between the two electron tracks, a, c, e) Phase 1 data; b, d, f) Phase 2 data. The 
statistical uncertainties on the data points are shown with error bars. 
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Figure 6.25: a, b) electron energy for external e'j events; c, d) photon energy and 
e, f) the cosine of the angle between the electron and the photon, a, c, e) Phase 1 
data; b, d, f) Phase 2 data. The statistical uncertainties on the data points are 
shown with error bars. 
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Internal background 


Channel 


Activity/mass (mBq/kg) 


208 fjj 


67 and 677 


10.1 ± 1.0 


212Bj 


from ^os^jj activity 


28.0 ±2.8 




from ^os^jj activity 


28.0 ±2.8 


207Bj 


67 


226.8 ±5.5 






51.3 ±2.4 




67 


24.6 ± 1.1 


234m 


l6 


47.7 ±3.0 




le 


213 ± 12 


214Bi 


6Ct 


3.35 ± 0.79 [58] 



Table 6.5: The internal backgrounds activities found from NEMO 3 data. 



To select external 67 events, the same selections are applied as described in Sec- 
tion 6.1.1 except that the external TOF hypothesis for an 67 event has to be greater 
than 4% and the internal TOF hypothesis has to be less than 1%. Figure 6.25 shows 
the electron energy, photon energy and cosine of the angle between the electron and 
the photon for Phase 1 and Phase 2 data. The number of external 67 events pre- 
dicted by simulation are 2480 ± 55 (stat) and 2399 ± 55 (stat), respectively. The 
number of data events are 2065 and 2016. The difference between the number of 
expected background events and data events gives 16.7% uncertainty on the external 
background in this channel. 



6.4 Summary and discussion 

The 67, 677, l6 and ea [58] channels were studied in order to measure the internal 
background activities with NEMO 3 data. Table 6.5 summarises the internal back- 
grounds, the channels studied and the activities measured. The activities of ^°^Bi 
and measured using NEMO 3 data are inconsistent with the HPGe detector 
result. The ^''^Bi activity is found to be 1.8 times more than HPGe measurement. 
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As is observed in Figure 6.5, the foil has regions with higher contamination of ^°^Bi 
than others and it is possible that HPGe measured a sample of ^^°Nd foil with lower 
contamination [83]. The activity of is four times higher than the HPGe mea- 
surement. The HPGe detector is a photon detector, whereas predominantly 
decays to an electron. The difference can therefore be due to the HPGe detector's 
inability to detect electrons. 

In this thesis a 17% uncertainty on the external background model is used from the 
validation of the external background model with the 67 channel. As the number of 
events from external background sources that mimic double beta decay is expected to 
be low, due to the small surface area of ^^''Nd foil (see Table 7.4), this uncertainty 
will not contribute significantly to the uncertainty on the measurement of 2i^/5/5 
decay half- life. 
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Chapter 7 



Two-neutrino double beta decay of 

This chapter describes the measurement of the 2z//5/5 half-hfe of ^^°Nd. The analysis 
is performed on data described in Section 5.4. The signature inside the NEMO 3 
detector for 2uj3j3 of ^^°Nd is two electrons emanating from this foil. The event 
selection criteria for two-electron events are given in Section 7.2 followed by the 
presentation of the half-life result in Section 7.3. 

7.1 Half- life definition 

The number of expected events, derived from the radioactive decay law and the 
event selection efficiency, e, of the considered channel, can be expressed as 

iV(t)=6iV,i(l-e"'°'^^\ (7.1) 

where Nat is the number of atoms in a sample. For the ^^''Nd sample in NEMO 3 
with a mass of 36.55 g, Nat is equal to 1.462 x 10^^. The time t is the data taking 
time, and T1/2 is the half-life of the studied decay mode. Because of the large value 
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Data 

2vpp 



a) 




b) 




Figure 7.1: The vertex distribution in the x — y plane (Ai? ) and the vertical coordi- 



nate (IA2I) for two electron events before the selections |Ai?| = \J (Ax)^ + (Ay)^ < 
2 cm and \^z\ < 4 cm are applied. The statistical uncertainties on the data points 
are shown with error bars. The 2vj3(3 signal MC (red) is scaled to the number of 
background subtracted data events. 



of Ti/2 for 2vf3(5 decay compared to t, Equation 7.1 can be rewritten as 



Ti/2 = eiVat In 2 



t 



N{t) 



eNat In 2 



t 



AT _ ATtot ' 



(7.2) 



where N^xp is the number of expected 2z//3/3 events found from data and Nl°l is the 
total number of background events remaining in the data after applying the event 
selection criteria. 



7.2 Two-electron event selection 

The following selections are applied to the data in order to reduce background to 
double beta decay: 

• Two negatively charged particles are required. The tracks must be associated 
with two scintillator hits. 

• Both tracks are required to originate from the ^^°Nd foil. 
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Figure 7.2: The distribution of the a) internal TOF hypothesis and b) external 
TOF hypothesis before applying the requirement of Pixlnt) > 0-04 and P{xlxt) < 
0.01. The external background events (gray) are normalised to the number of data 
events for -P(xLt) < 0-04 and P{xlxt) > 0.01. The error bars show the statistical 
uncertainties on the data points. 

• To remove the events from hot-spot regions, events are rejected if they originate 
from ^^''Bi and ^'^^™Pa contaminated regions defined in Equations 6.2 and 6.22. 

• The energy of each electron is required to be greater than 0.2 MeV. 

• Both electrons originate from a common vertex and have: 

- |Ai?| = y/{Axy + (Ay)2 < 2 cm, 

— \Az\ < 4 cm, 

where AR and Az are the distances between the two tracks' intersection points 
with the foil in the horizontal plane, x—y, and in the z coordinate. Figures 7.1a 
and Figures 7.1b show the distribution of these variables before applying the 
vertex cuts and after applying all other cuts described above. 

• In order to reduce the external background, the internal TOF hypothesis for 
two-electron events is required to be greater than 4%, and the external TOF 
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hypothesis to be less than 1%. Figure 7.2 shows the distributions of the in- 
ternal and external TOF hypotheses before applying the TOF hypothesis re- 
quirement and after applying all other selections described above. In these 
figures the number of external background events is normalised to the number 
of data events for P(xL < 0-04) and P(xL > 0-01). 

• The track length for each electron is greater than 30 cm. This cut value is 
lower than the one applied to channels involving one electron track only, as 
the TOF measurement for two track events is more precise [50]. Figure 7.3 
shows the track length distribution for one of the electrons before applying the 
track length requirement and after applying all other selections. It is observed 
that signal over background ratio is small for electrons with track length less 
than 30 cm and that the rate of events with track length <30 cm is not well 
simulated by the MC. 

• The tracks pass through one of the first two layers of the tracker. 

Table 7.1 gives the number of data events remaining after each event selection criteria 
is applied and Table 7.2 gives the fraction of MC events remaining after each cut. 
This shows the effect of cuts on the 2i^/3/3 signal and the internal and external 
backgrounds with major contributions in the 2z//5/5 signal region. 

In order to illustrate the effect of each cut on external backgrounds and radon in 
the tracker, the number of remaining simulated events after each cut is shown for 
^^^Bi in the glass of the PMTs and in the tracker. Figure 7.4 shows the display of a 
selected two-electron event originating from the ^^°Nd foil. The electron tracks and 
scintillator hits are shown. 
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20 40 60 80 100 1 20 140 160 180 200 



Electron track length[cm] 

Figure 7.3: The track length distribution of one of the electrons of two-electron 
events before applying the track length > 30 cm requirement. The error bars show 
the statistical uncertainties on the data points. 



NEM03 - RUN: 5142 EVENT: 41400 (2006-90-04 @ 22:07:32) 




Figure 7.4: A display of a two-electron event originating from ^^°Nd foil. The top 
{x — y) and side views {z coordinate) of the event are shown. 
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Cut 


Cut 


Number of events (N) 


(1) 


two tracks 


157575017 


(2) 


tracks associated with scintillators 


14309481 


(3) 


particles have negative charge 


1213245 


(4) 


tracks originated from ^^"Nd foil 


9249 


(5) 


Not originated from hot-spot regions 


8975 


(6) 


Energy of each electron > 0.2 MeV 


5657 


(7) 


AR < 2.0 cm and Az < 4.0 cm 


3983 


(8) 


Pixlt) > 0.04 and P(xLt) < 0.01 


2891 


(9) 


Length of each track > 30 cm 


2852 


(10) 


Tracks pass through one of the first two layers 


2789 



Table 7.1: Details of the two-electron event selection criteria. The number of events 
remaining in the data sample are presented after applying each selection. 



Cut 


2v(3P 




40k 


208 rp J 


228 Ac 


207Bi 


^^^Bi tracker 


2i4Bi PMT 


(1) 


31.0% 


0.6% 


1.9% 


0.55% 


0.30% 


0.07% 


0.0027% 


(3.5 X 10-4)% 


(2) 


15.6% 


0.3% 


0.08% 


0.27% 


0.15% 


0.035% 


0.0013% 


(1.8 X 10-4)% 


(3) 


13.6% 


0.19% 


0.036% 


0.21% 


0.10% 


0.026% 


0.001% 


(1.7 X 10-4)% 


(4) 


10.4% 


0.14% 


0.03% 


0.17% 


0.083% 


0.021% 


(3.37 X 10-4)% 


(3.9 X 10-5)% 


(5) 


10.1% 


0.12% 


0.028% 


0.16% 


0.079% 


0.020% 


(3.30 X 10-4)% 


(3.8 X 10-5)% 


(6) 


8.7% 


0.1% 


0.012% 


0.13% 


0.057% 


0.018% 


(3.37 X 10-4)% 


(3.8 X 10-5)% 


(7) 


7.8% 


0.08% 


0.0083% 


0.11% 


0.049% 


0.016% 


(1.23 X 10-4)% 


(2.5 X 10-5)% 


(8) 


7.4% 


0.078% 


0.0079% 


0.10% 


0.047% 


0.0152% 


(6.3 X 10-5)% 


(1.1 X 10-6)% 


(9) 


7.3% 


0.075% 


0.0078% 


0.10% 


0.046% 


0.015% 


(5.6 X 10-5)% 


(1 X 10-6)% 


(10) 


7.2% 


0.074% 


0.0077% 


0.10% 


0.0457% 


0.015% 


(4.9 X 10-5)% 


(1 X 10-6)% 



Table 7.2: The effect of cuts on 2i//9/9 signal and internal and external background 
MC samples. The fraction of events remaining after each cut is given for the internal 
backgrounds with major contribution to the 2v(3(3 signal and for ^^^Bi from the 
tracker and the PMTs. 
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7.3 Double beta decay half-life 



After applying all two-electron event selections, 2789 events remain. The number of 
Monte Carlo events are normalised to the activities found in Chapter 6 to estimate 
background remaining in the data sample. Due to the resolution of the tracking 
detector, events from neighbouring foils (^''^Mo, ^^Ca and ^^Zr) can be reconstructed 
as originating from the ^^'^Nd foil. The ^^Ca foil is heavily contaminated with 
which itself is a source of background to 2^(3(3 [84]. The activities of these isotopes 
are given in Table 7.3. The total number of background events from neighbouring 
foils is estimated to be 168.1 ± 13 (stat) ± 14 (syst) from MC simulations, where 
the systematic uncertainty is determined by the uncertainty on the half-life of these 
isotopes. 

Table 7.3 summarises the internal background originating from ^^'^Nd and neighbour- 
ing foils. The number of total internal background events is found to be 693 ± 28. 
The uncertainty on the number of background events is due to the statistical er- 
rors and the systematic uncertainties on the activity measurements. The systematic 
uncertainties on ^^'^Tl, ^^^Ac and ^^^Bi activities are fully correlated. The same is 
true for ^^^Pb and ^^^Bi. Table 7.4 summarises the external background that mimic 
2v(3(3 events. The total number of external background events is estimated to be 
53 ± 11. The 20% uncertainty includes statistical uncertainties and a 17% system- 
atic uncertainty on the external background model. By adding up the number of 
internal and external background events, the total number of background events is 
found to be 746 ± 30. 

The energy distribution of electrons with the lower energy, the energy sum of the 
two electrons and the opening angle between them are shown in Figure 7.5. The 
normalisation of the 2vf3l3 signal is found by scaling the simulated MC events to the 
same number of data events minus the radioactive background events. The data are 
in good agreement with the sum of the background and the 2^(3(3 signal distribution. 
The efficiency of the 2^(3(3 event selection is 7.2%. After background subtraction 
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Contaminant 


A (mBq) 


Selection efRciency 


Nbgr 


234mpg^ 


2.42 ±0.14 


0.074% 


144 ± 13 (stat) ± 9 (syst) 


40k 


10.8 ±0.6 


0.0077% 


66 ± 4 (stat) ± 4 (syst) 


208 J 


0.51 ±0.05 


0.10% 


46 ± 1 (stat) ± 4 (syst) 


228 Ac 


14.2 ± 1.4 


0.046% 


52 ± 2 (stat) ± 5 (syst) 


212Bi 


14.2 ± 1.4 


0.029% 


32 ± 2 (stat) ± 3 (syst) 


207Bi 


11.5 ±0.3 


0.015% 


138 ± 3 (stat) ± 5 (syst) 


214Bi 


0.17 ±0.04 [58] 


0.097% 


13 ± 1 (stat) ± 3 (syst) 


214pb 


0.17 ±0.04 [58] 


0.042% 


6 ± 1 (syst) 




2.6 ±0.1 


0.0094% 


19 ± 1 (stat) ± 1 (syst) 




1.25 ±0.05 


0.0087% 


9 ± 1 (stat) ± 1 (syst) 


48Ca 


(4.3 ±0.4) X 10-2 [84] 


0.062% 


2.1 ±0.2 (syst) 


goy 


30.6 ±3.1 [84] 


9.1 X 10-4% 


22 ± 3 (stat) ± 3 (syst) 


962r 


(7.2 ± 0.7) X 10-2 [85] 


0.082% 


5 ± 1 (syst) 




0.13 ±0.01 [86] 


1.34 X 10-3% 


139.0 ± 13 (stat) ± 14 (syst) 


Sum 






693 ± 28 



Table 7.3: Summary of the internal background to the 2^(3(3 signal of ^^''Nd. Each 
background's activity, the two-electron selection efficiency and the number of events 
which contribute to two-electron events in 924.7 days are given. The statistical 
uncertainty due to the finite size of MC samples are added in quadrature to the 
errors found in Chapter 6. 



Background 


^bgr 


Total external ^^^Bi 
Total external ^°Co 
2i4Bi/2i4pb in the tracker 
^^"^Bi in the tracker 


10.0 ± 4 (stat) 
5 ± 3 (stat) 

31 ± 3 (stat) 
7 ± 1 (stat) 


Sum 


53± 11 



Table 7.4: Summary of the external background remaining in the two-electron sam- 
ple after the event selection criteria in 924.7 days of data taking. A systematic 
uncertainty of 17% is added to the statistical uncertainties. 



110 




r\ 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
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Figure 7.5: The distributions of a) energy of the electron with lower energy; b) 
energy sum of the two electrons and c) the cosine of the angle between the two 
electrons. The sum of the total radioactive background (gray) and the 2i>/3(3 signal 
(red) is in good agreement with data. The statistical uncertainties on data points 
are shown with error bars. 

the half-life of the 2z//5/? decay is measured to be: 

T^r, = (9.11^°:! (stat.) ± 0.62 (syst.)) x lO^^ y. (7.3) 



Systematic uncertainties 

The breakdown of the systematic uncertainties on the 2z//3/5 half-life is as follows: 
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EsumEMeV] 



Figure 7.6: The distribution of energy sum of the two electrons for 2z//3/5 before and 
after smearing each individual electron energy by 2%. 

• The systematic uncertainty on the background estimate, obtained as the quadratic 
sum of the uncertainties on the estimates for the internal and external back- 
ground, is 4%, which translates into an uncertainty on the 2j//3/3 half-life of 
1.5%. 

• An uncertainty of 3% is added due to the uncertainty on the position of the 
iso^d foil in the detector. This is estimated by comparing the half-life results 
before and after changing the position of the ^^°Nd foil (Section 5.5). 

• The TOF cut is changed from 4% to 1% in order to estimate the uncertainty 
on this cut. This yields 1% uncertainty on the half-life results. 

• There is a 5% uncertainty on the efficiency calculation due to the inaccuracy of 
the GEANT simulation and tracking program. This was found by measuring 
the activity of the calibration sources [54]. 

• There is 2% uncertainty in the energy calibration coefficient [54] since the laser 
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energy calibration was not used. The uncertainty on the shape of the 2u(3(3 
energy sum distribution is studied by smearing the individual electron energies 
for 21/(3(3 simulated events by 2%. In this case the smeared energy is defined 
as 

El = E,{l + 0.02g), (7.4) 

where Ei is the energy of the i*'* electron and g is a. random number taken from 
a Gaussian distribution with a mean of zero and width of one. This yields a 
systematic uncertainty of 1.0% on the 21/(3(3 half-life. Figure 7.6 shows the 
distribution of energy sum of the two electrons for simulated 21/(3(3 events 
before and after smearing the electron energies by 2%. 



7.4 Comparison with previous results 

The 21/(3(3 decay of ^^''Nd was first observed in 1993 by a TPC experiment con- 
structed by the Irvine group in California [87]. The discovery was confirmed by 
the ITEP Moscow experiment later that year [88]. Both of these groups devel- 
oped their experiments in order to measure the half-life of this process. The TPC 
experiment in ITEP Moscow used 51.5 g of Nd203 and presented a half-life of 
Ti/2 = 1.88l[j|^ (stat) ± 0.19 (syst) x 10^^ y [89] in 1995. This half-life result 
had a large statistical uncertainty due to the short running time of the experiment 
(53 days). 

The Irvine group's TPC experiment was situated in an underground valve house in 
a canyon wall at the Hoover Dam in the USA [90]. The TPC was a rectangular 
polycarbonate box with internal dimension of 88 x 88 x 23 cm'^. The double beta 
decay source was in the centre and served as the central drift-field cathode. Anode 
and cathode wires were set perpendicular to each other and located near the walls 
of the TPC. A 1200 G magnetic field was applied perpendicular to the source plane. 
The helix of the particles provided the measurement of the momentum and thus 
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Figure 7.7: Comparison of NEMO 3 2uj3l3 half-life result (the result presented in 
this thesis) with ITEP and Irvine TPC experiments. The error bars present the 
quadratic sum of the statistical and systematic uncertainties. 

the energy. A Nd203 source with a mass 19.0 g enriched to 91% of ^^°Nd was used. 
After running the experiment for 262 effective days, 476 double beta decay events 
were observed. The efficiency of the detector for detecting two-electron events was 
13%, resulting in a measured half-life of Q.75toH (stat) ± 0.68 (syst) x 10^^ y [90]. 
Figure 7.7 compares the 21/(3(3 half-life result presented in this thesis with the results 
of the ITEP Moscow and Irvine TPC experiments. This measurement is slightly 
more than two standard deviations higher than the Irvine group results and two 
standard deviations lower than the ITEP Moscow result. 
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Chapter 8 

Limits on different modes of 
neutrinoless double beta decay 

This Chapter presents searches for beyond the Standard Model double beta decay 
processes as described in Sections 2.5.3 to 2.5.4. It starts by introducing the signa- 
ture for Oz//5/3 and OuPPx^ix^) signals in the NEMO 3 detector and continues with 
setting limits on the half-life of each individual mode using the method described 
in Section 5.8. Limits are also set on neutrinoless double beta decay of ^^'^Nd to 
excited states of ^^°Sm. 

8.1 Signal event selection 

The same two-electron event selections as described in Section 7.2 are applied to 
simulated events from: 

• Oz//5/3 decay to the ground state of ^^'^Sm via the mass mechanism (O^g (m^)). 

• Oz//5/? decay to the ground state of ^^°Sm via right-handed currents (Ogg RC). 

• Oz//3/3 decay with emission of Majoron(s). 
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0uf3(3 


Majorons 




0+ RC 


2+ RC 




n=l 


n=2 


n=3 


n=7 


efficiency (%) 


19.0 


10.9 


5.9 


2.1 


14.4 


12.2 


10.3 


4.9 



Table 8.1: Two-electron event selection efficiency for different 0^(3(3 decay modes. 



• OuPP decay to the 2f excited state {2f RC). The mass mechanism of this 
decay is highly suppressed due to the nuclear angular momentum and only 
the right-handed current is enhanced [91]. 

• OuPP decay to the Of excited state (O^*" {m^))- The right-handed current of 
this decay is highly suppressed due to the nuclear angular momentum and only 
the mass mechanism is enhanced [92]. 

The ffist three decays have the same event topology as 2upp. The event topologies 
of the last two decays are two electrons with emission of a photon (2^ RC) or two 
photons (0^ {rrii,)). In this thesis it has not been attempted to reconstruct the 
additional photon (s). 

Table 8.1 present the event selection efficiency for each of these modes. To set 
limits on the half-life of these modes, the whole energy sum distributions of the 
two electrons are investigated and no energy cuts are applied. Figure 8.1 shows the 
energy sum (Esum) distribution of different simulated modes, normalised to arbitrary 
numbers. In these figures the data and the sum of the backgrounds are also shown. 



8.2 Systematic considerations 

In all the above cases, there are energy regions where the 21/ PP background overlaps 
with the OuPP signal. This is especially true for the Majoron modes three and seven 
which have the largest overlapping energy regions with the 2uPP background. The 
2vPP background's normalisation factor (half-life) in Section 7.3 was found by fitting 
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Figure 8.1: Energy sum distribution of electrons for a) mass mechanism, right- 
handed current, decays to 2^ and 0"^ excited states and b) Majoron mode one 
(n=l), mode two (n=2), mode three (n=3) and mode seven (n=7). Data and 2v(3l3 
MC are also shown. The statistical uncertainties on the data points are shown with 
error bars. 

it to data minus radioactive background, assuming that there was no new physics 
in data. Thus, in order to have an unbiased signal search, the 2v(3(3 statistical 
uncertainty on the 2v(3(3 normalisation, which is found by fitting 2vj3(3 background 
to data minus radioactive background MC, was "floated" (see Section 5.8.3). This 
makes the 2^(3(3 background normalisation factor a free parameter. 

The uncertainty on the energy measurement affects the shape of the energy sum 
distribution of the signal and is different in each bin. Figure 8.2 shows the en- 
ergy sum of the two electrons for different simulated neutrinoless double beta decay 
modes before and after smearing the energy of each electron by 2% (Section 7.3). 
To find the shape uncertainty on the energy, the ratio (S-N) /N is plotted, where S 
is the smeared MC distribution and N is the nominal (unsmeared) MC distribution. 
These ratios for different Qv(3(3 modes are shown in Figure 8.3. The histograms 
show some statistical fluctuations. By fitting the histograms to different functions 
the statistical fluctuations can be reduced. Multi-Gaussian fits are performed for 
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Oz//?/3 modes and second order polynomial fits are performed for Oz//5/5 with emission 
of Majoron(s). Figure 8.4 shows the resulting histograms after the fit. The inverted 
histograms (shown in red) are also constructed in order to estimate the negative 
uncertainty. Figure 8.5 shows the shape uncertainty on the 2v(3(3 background. 
The uncertainties that are considered for setting limits on {]vf3(3 decays and Qvfifi 
decays with the emission of Majoron(s) are presented in Table 8.2. The shape sys- 
tematic uncertainties on signal and 2^(3(3 background and the efficiency systematic 
uncertainties on signal and radioactive background are assumed to be fully corre- 
lated. The shape systematic was found to contribute to the limit results with less 
than 1%. 



Uncertainty type 


signal 


2v(3l3 background 


other backgrounds 


Efficiency 


5% 




5% 


Activities 






4.0% 


half-life 




2.7% (Float) 




Energy measurement 


Shape 


Shape 





Table 8.2: The correlated systematic uncertainties on the ^^[3(3 signal, the 21/(3(3 
and radioactive backgrounds as considered for the limit setting. 

8.3 Limit Results 

8.3.1 Oiy(3f3 mass mechanism 

In order to search for OuPP, the Esum distributions of the Ou^P signal, 2z//3/3 and 
radioactive backgrounds are used as inputs to the limit calculation. The search is 
performed using the profile likelihood ratio test statistic described in Section 5.8.2 
(Equation 5.23). Table 8.3 gives the observed log likelihood ratio (LLRobs), the 
median expected log likelihood ratio (LLRmed) and the log likelihood ratio for the 
background with one and two standard deviations (LLRmcdilo" and LLRmediScr). 
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Figure 8.2: The distribution of energy sum of the two electrons for a) Ou/Sf] mass 
mechanism, b) Oz//5/5 right-handed current, c) OujSf] decay to 2^ excited state, d) 
OuPP decay to 0^ excited state, e) Majoron mode one, f) Majoron mode two, g) Ma- 
joron mode three and h) Majoron mode seven. The black lines show the energy 
distributions before smearing the energy of each electron and the red lines shows 
the energy distribution of the smeared MCs. 
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a) 
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d) 
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Figure 8.3: The distributions of (S-N)/N and fitted functions for a) Ou^f] mass 
mechanism, b). OuPP right-handed current, c) Ou^P decay to 2+ excited state, d) 
OuPP decay to O"*" excited state, e) Majoron mode one, f) Majoron mode two, g) 
Majoron mode three and h) Majoron mode seven. The error bars show the statistical 
uncertainties. 
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Figure 8.4: The estimated positive and negative shape systematics for energy distri- 
butions of a) Oz//5/3 mass mechanism, b) Oz//5/3 right-handed current, c) Oz//5/3 decay 
to 2+ excited state, d) OuPP decay to 0"*" excited state, e) Majoron mode one, f) 
Majoron mode two, g) Majoron mode three and h) Majoron mode 7. 
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Figure 8.5: The distributions of a) (S-N)/N and fitted second order polynomial 
function and b) the resulting positive and negative shape systematics for the energy 
distributions of 21/ f] (3 decays. The error bars show the statistical uncertainties. 



LLRobs 


LLRmed— 2cr 




LLRmed 


LLRmed + lo" 
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Table 8.3: Values of the observed and expected log likelihood ratios. 

The results in this table show that LLRobs = 4.6 is in agreement with LLRmcd within 
about one standard deviation and therefore there is no sign of a 0^(3(3 signal. The 
CLb value is 0.14 and therefore also consistent with a background only observation. 
To set limit on the observed and median expected number of events at 90% CL, the 
Oupp MC is re-scaled until CLs+b = 0.014 {CLs = 0.1). 

The observed upper limit on the number of events from the Ou^P mass mecha- 
nism is found to be 2.7 at 90% CL. With the detector efficiency of 19% and using 
Equation 7.2, the lower bound on the OuPP half-life is found to be: 



T°;2 > 1-8 X 10^2 y (90% CL). 



This is consistent with the median expected limit at 90% CL of 



T{72 > 1.33 X 10'' y (90% CL). 



12) 



This limit on the half-life is converted into a limit on the effective Majorana neutrino 
mass, (m^), using an NME of 3.14-4.04 [93] and G"'^ of 2.69 x 10-^^ y^^ [39]. The 
experimental lower limit on the half-life of ^^°Nd translates into an effective neutrino 



122 



mass limit of 

(m^) < 1.5 - 2.5 eV. (8.3) 

The uncertainty on the {niy) hmit is due to the uncertainty in NME calculations. 
Taking into account the nuclear deformation will modify this conclusion. The sup- 
pression of the NME for i^°Nd is estimated to be a factor 2.7 [94]. This mcreases 
the upper limit to 

(m^) < 4.0 - 6.3 eV. (8.4) 

Figures 8.6 shows the Egum distribution for Egum >2.5 MeV. The radioactive back- 
grounds, shown in grey, consist of 3.5 ± 0.9 ^os^l events and 0.6 ± 0.2 ^^^Bi events. 
Other backgrounds are found to be negligible in this energy region. A {)v(3(3 distri- 
bution with half-life of 1.8 x 10^^ y, corresponding to 2.7 events expected for the 
observed 90% CL, is shown in blue. 

8.3.2 Other Oi//3^ modes 



Table 8.4 gives the number of events corresponding to the median expected and 
observed limits and the corresponding lower bounds on the half-life of different 
neutrinoless double beta decay processes. All limits are set at 90% CL. Figures 8.7 
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Table 8.4: 90% median expected and observed upper limits on the number of events 
and lower limits on the half-life, T1/2 for different modes of neutrinoless double beta 
decay. 



and 8.8 show the Egum distribution of different neutrinoless double beta decay. Each 
signal is normalised to the same number of events as the corresponding exclusion 
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Figure 8.6: Distribution of the energy sum of the two electrons, Egum^ for Esum > 
2.5 MeV. The data are compared to the total background, consisting of radioac- 
tive background and the 2z//5/3 expectations. The error bars show the statistical 
uncertainties on the data points. A MC simulation of the signal with half-life of 
1.8 X 10^^ y is also shown. 

limit. As the 21/(3(3 background was left unconstrained, a different 21/(3(3 background 
normalisation is found for each exclusion signal limit. The lower limit on the half-life 
of the Majoron mode one leads to the upper limit on the Majoron coupling with the 
neutrino to be: 

{gMi) < (0.64 - 1.05) X 10"^ (8.5) 

To calculate this limit, the same NME as for the 0^(3(3 mass mechanism and the 
phase space factor, = 6.40 x 10~^^ [39] is used. Considering ^^°Nd defor- 
mation changes the limit to: 



{qmi) < (1.7-3.0) X 10- 



(8.6) 
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Figure 8.7: The Egum distribution of neutrinoless double beta decay to a) 2f and 
b) of excited states. The signal is normalised to the exclusion limit. The statistical 
uncertainties on the data points are shown with error bars. 
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Figure 8.8: The Egurn distribution of a) Majoron mode one, b) Majoron mode two, 
c) Majoron mode three and d) Majoron mode seven. The signal is normahsed to 
the exclusion limit. The statistical uncertainties on the data points are shown with 
error bars. 

The upper limit on the coupling of neutrino with right-handed currents and Majoron 
mode two and three are not presented in this thesis as the NME calculations of these 
modes are out of date [95]. 
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8.4 Summary and discussion 



The observed limit on the half-hfe of neutrinoless double beta decay of ^^''Nd was 
found to be 

r°;2 > 1-8 X 10^^ y (90%cl). 

This half-life limit has improved the previous limit on this isotope by a factor of ten. 

The best previous limit results were measured by the Institute for Nuclear Research 
of the USSR in 1986 [96] for 50.5 g of enriched ^^''Nd. The experiment had a sim- 
ple set-up using a central foil surrounded by four scintillation counters providing 
calorimetry and event reconstruction. This experiment gave a lower bound on the 
Oz//3/9 half-hfe of, T^^^ > 1.7 x 10^^ y (95% CL). Despite this significant improvement 
on the half-life limit, the upper bound on effective neutrino mass has a large uncer- 
tainty due to the uncertainties on the NME calculations compared to other Ou/SP 
isotopes. Therefore, further progress in the calculation of the NME for Ou/Sf] decay 
of ^^°Nd is required to improve the {m,^) limit found from this isotope. The NME of 
other neutrinoless double beta decay modes are also required to be recalculated with 
improved calculating tools. This is especially true for right-handed current Ou^P 
searches. 

The lower limit on the half-life of the Majoron mode one decay of ^^°Nd was measured 
to be 

T°f > 1.6 X 10^^ y (90%CL), 

which has improved the previous limit significantly (T°/f > S.SlxlO^o y (90%CL) [90]). 
The corresponding limit on the Majoron-neutrino coupling is 

{gMi) < (1.7-3.0) X 10-^ 

This limit is comparable with limits found with ^''^Mo {{gMi) < (0.4 — 1.9) x 10~^)) 
and ^^Se ((0.66 — 1.7) x 10~^) [55]. This result is obtained although only 36.55 g of 
^^°Nd is used, in comparison with 6.9 kg of ^""Mo and 0.93 kg ^^Se. This is due to 
the dependence of on Q^^. 
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Chapter 9 



Conclusion 



In this thesis a precise measurement of the 2i^/3/3 half-hfe of ^^°Nd and searches for 
different modes of neutrinoless double beta decay using NEMO 3 data have been 
presented. These involve measurements of the activities of backgrounds originating 
from the internal sources and validation of the external background model. The 
results given in this thesis have been submitted for publication [97]. 

The internal background activities are measured by studying 67, 677 and le chan- 
nels. Hot-spot regions due to ^'^^Bi and ^^^'"Pa contamination are observed in the 
67 and le channels, respectively. These regions are subsequently removed from the 
analysis. The ^°^Bi and activities found in this thesis have invalidated the HPGe 
results for these isotopes. The ^'^^™'Pa, ^^''Bi, ^°K, ^^^Ac isotopes are found to be the 
major internal background sources to 2i'(3f3 and their activities are measured with 
uncertainties less than 6%. The ^os^l isotope is the major background to Oz//3/3 and 
its activity is measured with an uncertainty of 10%. 

The number of two-electron background events originating from the neighbouring 
double beta decay sources and their contaminants are found to be 22% of the total 
background events. The external background model used in NEMO 3 was validated 
by selecting external particles that are interacting with ^^'^Nd foil. A 17% uncer- 
tainty on the model is used from the validation with the external 67 channel. The 
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total external background contribution to 2^(3(3 is found to be 7% of the total back- 
ground events. After applying the two-electron event selection and subtracting the 
background events, the half-life of 2^(3(3 is measured to be: 

Tlj^ = (9.1ll°:l (stat.) ± 0.62 (syst.)) x 10^« y. 

The dominant systematic uncertainties are due to the uncertainties on the tracking 
efficiency (5%), ^^°Nd foil position (3%) and number of radioactive background 
events (1.5%). 

A search is performed for different Oz//3/3 modes by studying the energy sum distri- 
butions of the two electrons in the final state. Limits on Qvf3l3 signals are set using 
the profile-likelihood ratio technique and by floating the 2^(3(3 background. This 
limit setting technique has been used for the ffist time in the NEMO experiment. 
The lower limits on the neutrinoless double beta decay half-lives of ^^°Nd for the 
mass mechanism (Oz//5/?) and the Majoron mode one decay {Qv(3(3x) are, 

T°;2 > 1-8 X 10^^ y (90%CL), 
T^/f^^ > 1.6 X 10^1 y (90%CL). 

Both results significantly improve the previous limits on neutrinoless double beta 
decay half-lives of this isotope and lead to the upper limits on the effective neutrino 
mass and Majoron- neutrino coupling of 

(m^) < 1.5- 6.3 eV, 
{gMi) < (0.64 - 3.0) X 10"^ 

The NME calculation of Oz//?/3 must be improved in order to reduce the uncertainty 
on these limits. Limits are also set for the first time on half-lives of right-handed 
currents, excited state Oz//3/3 decays and several other models leading to emissions 
of Majoron(s). 
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